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ABSTRACT 

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data 
from the Sloan Digital Sky Survey-Ill: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). 
From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i < 21.8 quasars 
over an area of 2236 deg 2 with confirmed spectroscopic redshifts between 2.2 < z < 3.5, filling in 
a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness 
of the survey through simulated quasar photometry, and check this completeness estimate using a 
sample of quasars selected by their photometric variability within the BOSS footprint. We investigate 
the level of systematics associated with our quasar sample using the simulations, in the process 
generating color-redshift relations and a new quasar ^-correction. We probe the faint end of the QLF 
to Mi(z = 2.2) ?» —24.5 and see a clear break in the QLF at all redshifts up to z — 3.5. We find that a 
log- linear relation (in log$* — M*) for a luminosity and density evolution (LEDE) model adequately 
describes our data within the range 2.2 < z < 3.5; across this interval the break luminosity increases 
by a factor of ~2.3 while $* declines by a factor of ~6. At z < 2.2 our data is reasonably well fit 
by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", 
in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to 
a number of theoretical models and find that models making a variety of assumptions about quasar 
triggering and halo occupation can fit our data over a wide range of redshifts and luminosities. 
Subject headings: surveys - quasars: demographics - luminosity function: AGN evolution 
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1. INTRODUCTION 

Quasars, i.e. luminous active galactic nuclei (AGN), 
represent a fascinating and unique population of objects 
at the intersection of cosmology and astrophysics. The 
cosmological evolution of the quasar luminosity function 
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(QLF) has been of interest since quasars were first iden- 
tified a half-century ago (lSandagdll96lt Hazard et al 
l96l ISchmidtl 11963: [Qkill963t IGreenstein fc Matthews 



19631: lBurbidgeHl967l ). 

Measuring the QLF, and its evolution with red- 
shift, is important for several reasons. It is gener- 
ally believed that present-day supermassive black holes 
(SMBHs) gained most of their mass via gas accre- 
tion during an active nuclear phase, potentially at 
quasar luminosities (Xp n i > 10 45 erg s -1 



| Salpeterl 11961: 

Zel'dovich fc Novikovi 119651 : iLvnden-Belf 119691 : iSoltanl 



1982T ) . so an accurate description of the QLF allows us 



to place constraints on the formation history of super- 
mass iv e black holes ( e. g.. iReesl Tl984t iMadau fe Reesl 
20011 iVolonteri et all 120031: iVolonteri fc Reesl 120061: 



Netzer fc Trakhtenbrotll2007l : lHaimadl2012l ) and to map 



the black hole accretion history of the Universe v i a the 
black hole mass fun ction dShankar et al.ll2009l [20Tot ISher] 
120091: IShen fc Kellvl[2012]) . as well as con strain the effect 
of black hole spin on the ce ntral engine (|Volonteri et al.l 
[20051: IFanidakis et al.l[20H . 

Measurements of the QLF also place constraints on 
the intensities and nature of various cosmic back- 
grounds, including t he buildup of the c o smic X-ray 
(iShanks et al.l 119911: iComastri et~al1 119951: lUeda et al. 
2003[ iBrandt fc Hasingerl 120051: iHickox fc Markevitch 
20061). ultraviolet (iHenrvl 119911) and infrared (IR) 
(jHauser fc Dwekl 120011 : iDole et all 120061 ) backgrounds. 
Knowledge of the UV background is relevant for calcu- 
lations that involve the contribution of quasar UV pho- 
tons to the epoch of H reionization (e.g-. lFan et al.|[2006l) 
at z > 6. At lower (z < 6) rcdshift, quasars con- 
tribute towards a fraction of the ionizing photons that 
keep most of the H ionized, allowing studies of the Ly-a 
forest (LyaF; e.g. lLvndslll97ll : lMeiksinll2009| ). Helium 
reionization ( He II— »He HI) can be measured by its effect 
on the LyaF jlJakobsen et al.lll994t iReimers et allll997l: 
Smette et all 120021: IReimers et al.l 120051: iSvphers et al.1 



20111: iWorseck et alj|2011l) . This second epoch of reion- 
ization occurs at z ~ 3, and may be driven by UV pho- 
tons from quasars, so an accurate determination of the 
QLF at this epoch is a key consistency check on the He 
reionization measurements. 

Furthermore, the co-evolution of galaxies and AGN is 
a crucial ingredient in, and test of, modern theories of 
galaxy formation. The energy feedback from AGN is 
thought to impact their host galaxies, an d thus influenc- 
ing their present-d ay properties (e.g., iCattaneo et al.1 
120091: iFabianl I2012D . Observations of the evolution of 
quasar properties over cosmic time can inform such mod- 
els and therefore our understanding of the galaxy-black 
hole connection. 

Recent large quasar surveys have allowed us to 
study the properties of the quasar population with 
unprecedented statistical precision. The number of 
known quasars has increased nearly 100-fold since the 
late 1990s, (for photom etrically identified quasars, 
see iRichards et al.1 120091 ) and since that time, there 
has been a larg e effort to measure the QLF in 
the UV /o p tical dBovle et al 
Wolf et all JgfjM IHunt et al 



2001: IJiangeTaLl I2008L [2001 ICroom et al.l [2 009a: 



Glikman et all 120101: iWillott et al l 120101; iGlikman et all 



20111: llkedaet all I2011L I2012h IMasters et all I 2~ 
mid-infrare d dBrown et all 120061 : ISiana et"all 120081 : 
lAssef et all 1201 ll) and the s o ft and hard X-ray 
dCowie et all 120031: lUeda et al.l 120031: lHasinger et all 
20051: iBarger et all [20051: iSilverman et al.l I2005L [ 20081: 
Aird et all 20081; iTreister et al.l 120091: lAird et all bold: 
Fiore et al. 2013). An overview of recent determinations 



of the optical QLF is given in Table [TJ 

Qu asar number density evolve s strongly with red- 
shift (ISchmidtJ 119701: lQ~3menfl98"l [Schmidt et all[l995l: 
Fan et al.l 120011: IRichards et all l2006bl: ICroom et all 
2009bfl . and one of the key goals of quasar studies is to 
understand what drives this strong evolution. A caveat 
here is that the evolution of the optical QLF is a com- 
posite of intrinsic quasar evolution and the evolution of 
the obscuring medium in quasar hosts. In this study, we 
concentrate on the unobscured quasar population, de- 
fined as objects that were selected via their UV/optical 
rest-frame continuum and the presence of broad, >1000 
km s — , emission lines. We leave investigations of 
the obscured AGN po pulation to other studies, e.g., in 
the m i d-infrared (e.g.. lLacv et al.l 120041: IRichards et all 
I2006at [Stern et al.l l2012t lYan et alll2012D. and X-ra 



(iTueller et alll201(|lLuo et all 1201 lHCorral et all l2(Hl[ 
Brig htman fc Uedai [20121: iLehmer et all 1201 2l)" ~~As the 

QLF is observed to have a broken power-law form, it is 
necessary to probe below the luminosity at which the 
power-law breaks in order to distinguish luminosity evo- 
lution (where the luminosity of AGN changes with time, 
but their number density remains constant) from density 
evolution (where the number density of AGN changes, 
but the luminosities of individual objects remains con- 
stant), or a combination of the two. 

The QLF is defined as the number density of quasars 
per unit lumino sity. It is often d escribed by a dou- 
ble power-law dBovle et all 120001: ICroom et ail 120041: 
IRichards et alll2006bl hereafter. R06) of the form 



HL,z) = 



{L/L*) a + (L/L*Y 



(1) 



with a characteristic, or break, luminosity L*. An al- 
ternative definition of this form of the QLF gives the 
number density of quasars per unit magnitude, 



lQOA(a+l)[M-M"(z)] _j_ lQDA(p+l)[M-M*(z)] 



(2) 

The dimensions of $ differ in the two conventions. We 
have followed R06 such that a describes the faint end 
QLF slope, and /3 the bright end slope. The a/0 con- 
vention in some other works (e.g.. [Croom et a"Ll r2009a) 
is in the opposite sense from our definition. Evolution 
of the QLF can be encoded in the redshift dependence 
of the break luminosity, </>*, and also potentially in the 
e volution of the pow er-la w slopes. 

iBovle et al.l (pOOOT ) and ICroom eTall (pOOl found that 
the QL F measured in the 2dF Quasar Redshift Survey 
(2QZ, ICroom et all 120041) was well fit by a pure lumi- 
nosity evolution model where $i ^ was constant but 
Fontanot et al.ll2007l : iBongiorno et al.ll2007l : iReves et all M* evolved with redshift. In this model, M* changed 



Croom et all 120041: lHao et all 120051: IRichards et al l 



2000; 



2004; 



Fan et all 



Fan et all 



2001 



2004; 



20051 I2006bl: iFan et all l200a I Jiang et al 
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TABLE 1 

Selected optical quasar lumino sity function measurem ents. 
a Cosmic Evolution Survey iIScoville et al.||2007bI) . 
b No Type-1 quasars were identi fied, though a low-l uminosity z ~ 5.07 Type-2 qu asar was discovered. 
c NOAO Deep Wide-Field Survey [Jannuzi fc Dey 1999) and the Deep Lens Survey (IWittman et al.| [2002). 

d SDSS Faint Quasar Survey. 
e The "boss21" area on the SP SS Stripe 82 field. 
^2dF-SDSS LRG And QSO Survey (ICroom et alJI2009bB . 
9 Photometric sample from SDSS; spectroscopic confirmation from SDSS and other telescopes. 
h Canada-France High-z Quasar Survey (Willott et al. 2003) 
; 2dF Quasar Redshift Survey iICroom et AL.II20O4T) . 
3 From our "uniform" sample defined in Section 12.31 
fe from a catalog of >1, 000, 000 photometrically classified quasar candidates. 



Survey 


Area (deg^) 


Nq 


Magnitude Range 


z-range 


Reference 




OOODS(+SDSS) 


0.1+(4200) 


13(+656) 


22.25 < z 8 50 < 25.25 


3.5 < z < 5.2 


Fontanot et al. (2007) 




VVDS 


0.62 


130 


17.5 < 7 A B < 24.0 


< z < 5 


Bongiorno et al. (2007) 




COMBO-17 


0.8 


192 


R < 24 


1.2 < z < 4.8 


Wolf et al. (20031 




COSMOS" 


1.64 


8 


22 < i' < 24 


3.7 < z < 4.7 


Ikeda et al. 72011 ) 




COSMOS 


1.64 


b 


22 < i> < 24 


4-5 < z < 5.5 


Ikeda et al. (2012) 




COSMOS 


1.64 


155 


16 < I AB < 25 


3 < z < 5 


Masters et ah^2012) 




NDWFS+DFS C 


4 


24 


R < 24 


3.7 < z < 5.1 


Glikman et al. 75571) 




SFQS d 


4 


414 


g < 22.5 


z < 5 


Jiang et al. (2006) 




BOSS e +MMT 


14.5+3.92 


1 877 


9 $ 23 


0.7 < z < 4.0 


Palanquc-Dclabrouillc et al. 


(2012) 


2SLAQ/ 


105 


5 645 


18.00 < g < 21.85 


Z < 2.1 


Richards et al. (2005) 




SDSS 9 


182 


39 


i < 20 


3.6 < z < 5.0 


Fan et al. (2001) 




SDSS+2SLAQ 


192 


10 637 


18.00 < g < 21.85 


0.4 < z < 2.6 


Croom et al. (2009a) 




SDSS Main+Deep 


195 


6 


z AB < 21.80 


z ~ 6 


Jiang et aT"72009) 




BOSS Stripe 82 


220 


5 476 


i >18.0 and g <22.3 


2.2< z <3.5 


PalaiiQuc-Dclabrouillc et al. 


(2011) 


CFHQS' 1 


500 


19 


z' < 22.63 


5.74 < z < 6.42. 


Willott et al. (2010) 




2QZ* 


700 


23 338 


18.25 < 6j < 20.85 


0.4 < z < 2.1 


Bovle et al. (2000); Croom et al. (2004) 


SDSS DR3 


1622 


15 343 


i < 19.1 and i < 20.2 


0.3 < z < 5.0 


Richards et al. (2006b) 




BOSS DR9 


2236 


J 23 201 


g <22.00 or r <21.85 


2.2< z <3.5 


this paper 




SDSS DR7 


6248 


57 959 


i < 19.1 and i < 20.2 


0.3 < z < 5.0 


Shen & Kelly (2012) 




SDSS Type 2 


6293 


887 


Lota > 10 8 ' 3 L o 


z < 0.83 


Reyes et al. (2008) 




SDSS DR6 fc 


8417 


> 850,000 


t < 21.3 


z ~ 2 and z ~ 4.25 


Richards et al. (2009) 





from « —26.0 to m —22.0 between z ~ 2.5 and z ~ 
0. However, this paradigm is ch allenged using recent, 
deeper data. ICroom et al.l (|2009aT ) measured the optical 
quasar luminosity function at z < 2.6 from the combi- 
nation of the 2dF-SDSS LRG And QSO survey (2SLAQ; 
ICroom eT al. 2009rJ), which probes down to a magnitude 
limit of q = 21.85, and the SDSS-I/II Qu asar survey 
(|Richards et all [20021 : ISchneider etldl l27jTol ) to i = 19.1 
(z < 3) and i — 20.2 (z > 3). Here, the double power-law 
form with pure luminosity evolution provides a reason- 
able fit to the observed QLF from low z up to z ~ 2, but 
it appears to break down at higher redshift. However, the 
2SLAQ sample has few objects above z ~ 2, and SDSS 
does not probe down to at higher rcdshifts, making it 
difficult to constrain the faint end of the QLF at high z. 

At z > 2, the constraints on the QLF are less clear- 
cut, as the selection of luminous quasars becomes less 
efficient. This situation arises because the broad-band 
colors of z ~ 2.7 and z « 3.5 quasars are very similar 
to those of A and F stars (IPanl fl999t iFan et all 120011 : 
iRichards et "all [2001 IRoss et al.ll2012D in the Sloan Dig- 
ital Sky Survey color system {Fukugita etallH99^ . Al- 
though we have goo d constraining power at the bright 
end at z > 2, (e.g. IRichards et al.1 [2006bt Uiang et all 
2009), there is uncertainty in the form, and evolution of 
the QLF at z > 2, especially at the faint end. The red- 
shift range z ~ 2 — 3 is of particular importance since 
the luminous quasar number density peak s here; this is 
often referred to as the "quasar epoch" (|Osmerl 119821 : 



Warren et al1ll99l I Schmidt, et alJllflM IFan et al.1 12001: 
Richards et alJl2006bt ICroom et al.ll2009al ). 



For our study, we use data from the SD SS-III: Baryon 
Oscil lation Spectroscopic Survey (BOSS; iDawson et al.l 
2012) that is specifically designed to target fa int, q < 22, 
quasa rs in the redshift range z — 2.2 — 3.5 (Ro ss et alJ 
[2012] ). The first two phases of the SPSS C' SDSS- 
I/IF, hereafter simply SDSS: lYork et al.1 12000) have 
been completed (|Abazaiian et al.ll2009D . with a sample 
of » 10 0, 000 spectroscopically confirmed quasars at < 
z < 5 ([Schneider et al.l 12010 ). The thir d incarnation of 
the S loan Digital Sky Survey (SDSS-III; lEisenstein et al.1 
2011) is taking spectra of 150,000 z > 2.2 quasars as 
part of the BOSS. The main scientific motivation for the 
SDSS-III BOSS Quasar survey is to measure the baryon 
acoustic oscillation feature (BAO) in the Lyman-a forest 
(LyaF; iSlosar et al.l 1201 TT ) . This sample is designed to 
select quasars with 2.2 < z < 3.5, and will have an order 
of magnitude more objects at z > 2 than SDSS, sampling 
the quasar luminosity function ~ 2 magnitudes deeper at 
each redshift. Combining the BOSS and SDSS observa- 
tions gives a dynamic range of ~5 magnitudes at a given 
redshift, and a primary mo tivation for our s tudy is to ex- 
tend the work presented in [Richards et al.1 ([2006b| ). both 
in dynamic range in luminosity, and concentrating on 
the redshift range z = 2.2 — 3.5, where the original SDSS 
selection was sparse-sample d in an attempt to m inimize 
the contamination by stars (Ri chards et al.ll2~002j) . 
In this paper we present the optical quasar luminos- 
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ity function (QLF) from the first two years of BOSS 
spectros copy, data include d in SDSS Data Release Nine 
(DR9; lAhn et all 12012ft 27 . We use data from the 
3671 deg 2 observed over the DR9 footprint, and sup- 
plement this with deeper data over a smaller area (14.6 
deg 2 ), in order to probe the redshift range 0.7 < 
z < 2.2, also observed as part o f the BOSS (Table [TJ 
iPalanque-D clabro uille et al.l 120121) . Table [1] places the 
BOSS DR9 survey in context as a wide-field, medium- 
depth survey, and we will return to the surveys that 
match BOSS in redshift. 

The outline of this paper is as follows. In Section 2 
we describe our data sets, which includes both color and 
variability selected AGN samples. In Section 3, we quali- 
tatively compare our different quasar samples, and quan- 
tify our selection function using both empirical data, and 
new, updated template quasar spectra. In Section 4 we 
present the SDSS+BOSS quasar number counts and a 
new quasar fc-correction based on our simulations. In 
Section 5, we present the combined SDSS+BOSS QLF, 
sampling the range —24.5 < Mj < —30 in absolute mag- 
nitude across redshifts 0.7 < z < 3.5 and compare to 
previous measurements. In Section 6, various models of 
the double-power law form are fit to our data, we com- 
pare our results to recent theoretical predictions in the 
literature, and place our new results in a broader con- 
text. We present our conclusions in Section 7. In Ap- 
pendix we investigate further the selection function 
models introduced in Section 3 and in Appendix [B] pro- 
vide tables of our measured QLFs. For direct comparison 
with, and extension of, R06, we assume a flat cosmology 
with f2 A =0.70 and H = 70 h' 1 Mpc. O ur magnitudes 
are b ased on the AB zero-poin t system (lOke fc Gunnl 
119831 ) and are PSF magnitudes (iStoughton et al.1 120021)7 
correc ted for Galactic extinction following iSchlegel et al.1 
(1998). Absolute magnitudes (M) are det ermined using 
luminosity distances for this cosmo logy (Peebles 1980, 
[i99l IHogg eTaII[2002l: IWriglrdl2006D . 

2. DATA 

We use imagin g data that are p art of Data Re- 
lease Eight (DR8; lAihara et al.ll2011f) in order to select 
spectroscopic tar gets that form the D a ta Re lease Nine 
(|Ahn et al.l [2012ft dataset. iRoss et al.l (|2012D describes 
the BOSS quasar target selection algorithms used to 
identify objects for spectroscopy. In summary, we use the 
subset of the DR9 data that employs the "Extreme De- 
convolution" (XDQSO) algorithm of IBovv et all (|201lD 
to select quasars based on their optical fluxes and col- 
ors to define a uniform sample. The XDQSO procedure 
is supplemented by a selection using optical variability, 
where we have repeat imaging data within the DR8 foot- 
print. The DR9 data are the first two years of BOSS 
spectroscopic data, a nd the full DR9 Q uasar Catalog 
(DR9Q) is detailed in iParis et all (|2012l ). Fig Q] shows 
the sky coverage of the DR9 quasar dataset. However, 
the XDQSO selection was not implemented in the first 
year, leading to effects on completeness that we will ad- 
dress below in order to perform a QLF measurement. 

2.1. Imaging and Target Selection 
27 http://www.sdss3.org/surveys/ 




40° 20° 0° -20° -40° 

RA (degrees) 



Fig. 1.— The sky coverage of the SDSS-III: BOSS DR9 quasar 
dataset (colored regions) overlaid on the final expected footprint of 
BOSS (gray). These areas are 3671 and 10 269 deg 2 , respectively. 
The upper panel shows the coverage in the NGC, and the lower 
one is the SGC. Each sector (covered by a unique combination of 
spectroscopic tiles) is colored according to the fraction of quasar 
targets, selected with the uniform XDQSO algorithm, which have 
successful redshifts. The sectors which contribute data to the QLF 
analysis have > 85% spectroscopic completeness (yellow-orange- 
red regions). This area is 2236 deg 2 . The Stripe 82 field runs from 
-43° < R.A. < 45° at Decl. = ±1.2 5° and generally has / sc < 0.85 
(where / sc is defined in Sec. I3.l!4l . However, since Strip e 82 had 
quasar targets that were variability selected (see Sec. 12.411 . the true 
number of quasars in this field is actually very high. 



The SDSS-IILBOSS uses the imaging d ata gathered 
by a dedicated 2.5m wide-field telescope (jGunn et al.l 
2006), which collected light f or a camera with 30 
2kx2k CCDs dGunn et al.l I1998D over five broad bands 
- ugriz (FukugitaeLaU 11996) - in order to image 14,555 
unique deg 2 of the sky. This area includes 7,500 deg 2 
in the North Galactic Cap (NGC) and 3,100 deg 2 
in the South Galactic Cap (SGC). The imaging data 
are taken on dark p hotometric nights of good see- 
ing (IHogg et al.l I2001D and are calibrated photometri- 
cally ([Smith et al.l 120021: llvezic et al.l 120041: iTucker et all 
120061: iPad manabh an et al.l 120081 ). and astrometrically 
(iPier et al.ll20(H) before object param eters are measured 
jLupton et al.ll2001t IStoughton et al.ll2002ft . 

Using the imaging data, BOSS quasar target candi- 
dates are selected for spectroscopic observation based on 
their fluxes and colors in SDSS bands. However, selec- 
tion of quasars for BOSS spectroscopy is complicated by 
two facts: (i) The optic al colors of z ~ 2.7 quasars resem- 
ble faint A and F stars (|Fadll99ftlRichards et alj|2001bt 
IRoss et al.H2012ft and (ii) to maximize the number den- 
sity of quasars for LyaF cosmology, we are required to 
work close to the magnitude limit of the (single-epoch) 
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imaging data, leading to larger photometric errors, ex- 
pansion of the stellar locus and higher stellar contam- 
ination. All objects classified as point-like and having 
magnitudes of g < 22 or r < 21.85 are passed to the 
quasar target selection code. 

As was the case for the original SDSS Quasar survey, 
radio data was used to select quasars. Specifically, op- 
tical stellar objects with g < 22 or r < 21.85 which 
have matches within 1" to radio sources apparent in 
the Faint Radio So urces at Twenty cm (FIRST) survey 
(|Becker et al.l Il995l ) are considered as potential quasar 
targets, irrespective of their radio morphology. Approx- 
imately 2% of targe ts, a nd wl.3% of our uniform quasar 
sample (defined in §[2]3]below), satisfy the radio selection 
criteria. 

As the main science goal of the BOSS quasar sam- 
ple is to probe the foreground hydrogen in the IGM, 
priority was placed on maximizing the surface den- 
sity of z > 2 quasars (jMcDonald fc Eisensteinl 120071: 
iMcQuinn fc White] 1201 ID . rather than creating a homo- 
geneous dataset. The target selection is consequently a 
compl icated heterogen ous combination of several meth- 
ods (jRoss et aLll2012f). However, a uniform subsample 
(called "CORE" in iRoss et al.l l2012) was defined to al- 
low statistical studies of quasar demographics to be per- 
formed. The spectroscopic observations, and creation of 
this uniform subset of objects, are described in the next 
two sections. 



i i i i i i 




redshift 



Fig. 2. — Quasar N(z) redshift distributions. The dotted red 
histogram shows the redshift distribution for the full SDSS-III: 
BOSS DR9 quasar dataset, while the solid red line shows those 
objects uniformly selected by the "XDQSO" method across 2.2 < 
z < 3.5. The bl ack histogram is the fi nal distribution from the 
DR7Q catalog of Schneide r et al.l pOlOl '). 



2.2. Spectroscopy 

The BOSS spectrographs and their SPSS p redecessors 
are described in detail by iSmee et al.l ((2012) . In brief, 
there are two double-armed spectrographs that are sig- 
nificantly upgraded from those used by SDSS-I/II. They 
cover the wavelength range 3600 A to 10, 400 A w ith a 
resolving power of 1500 to 2600 (|Smee etal.1 12012ft . In 
addition, the throughputs have been increased with new 
CCDs, gratings, and improved optical elements, and the 
640-fibre cartridges with 3" apertures have been replaced 
with 1000-fibre cartridges with 2" apertures. Each ob- 
servation is performed in a series of 900-second expo- 
sures, integrating until a minimum signal-to-noise ratio 
is achieved at a fiducial magnitud e for the given spectro- 
scopic plate (|Dawson et al.ll2012l ). 

Once target selection is completed, the spectroscopic 
targets are assigned to tiles of diameter 3° using an al- 
gorithm that is adaptive to the density of targets on 
the sky (|Blanton et alj 120031 ). Of the 1000 available 
fibers on each tile, a maximum of 900 fibers are al- 
located for science targets, of which ^160-200 are al- 
located to quasar targets, while 560-630 fibers are as- 
signed t o galaxy targets, an d 20-90 to ancillary science 
targets ([Dawson et al] l2012). Because of the 62" diame- 
ter of the cladding around each optical fiber, two targets 
with a separation smaller than that angle cannot both 
be observed on a given spectroscopic plate, and differ- 
ent classes of targets are assigned priorities when such 
a collision arises. CORE quasars are assigned higher 
tiling priority th an the galaxy targets (Appendix B of 
IRoss et al.ll2012T) . To cover the survey footprint without 
leaving gaps, adjacent tiles overlap, alleviating the fiber 
collisions problem somewhat. This leads to the definition 
of a "sector" - a region covered by a unique set of tiles (see 



Blanton et"aT]|2003l; iTegmark et aLll2004i: ISwanson et al.l 
2008HWhite et al.ll201lD . As in previous SDSS analyses, 
we work on a sector-by-sector basis to define our various 
completenessess. 

The DR9 footprint is 3 671 deg 2 , and is given in Fig.[TJ 
In total, we obtained 182 973 spectra of objects that were 
selected as BOSS quasar targets, and lBolton et al.l (|2012j ) 
describes the automated spectral classification, redshift 
determination, and parameter measurement pipeline 
used for the BOSS. A total of 167 331 of these had the 
specPrimary flag set to 1, indicating that this was the 
best spectroscopic observation of an object; this cut, 
by definition, re moves objects with duplic a te sp ectra. 
As described in lAdelman- McCarthy et ail (|2008l ) and 
iBolton et all (|2012f ). each redshift is accompanied by a 
flag, zWarning, which is set when the automatically de- 
rived (a.k.a. pipeline) redshift and classification are not 
reliable; 132 290 of these objects do not have this flag set. 
Of these, 54 019 have pipeline redshifts between 2.2 and 
3.5. A summary of these numbers is given in Table [2] 

Each of the quasar target spectra has also been vi- 
sually inspected, and the redshift corrected where nec- 
essary. In total there are 87 822 spectroscopically con- 
firmed quasars in the DR9Q, whi le approximately half 
the quasar candidates were stars (| Paris etaDEOli). If 
there was confidence in a secure redshift from the visual 
inspection of the spectrum, the flag z_conf .person was 
set to be > 3. Most of the zWarning ^ objects have 
secure redshifts after visual inspection; over 97% of the 
specPrimary BOSS quasar target spectra have secure 
redshifts (Table [5]). Among these objects, 54 593 are 
confirmed, by visual inspecti on, to be at 2.2 < z < 3.5. 
For comparison the DR7Q ([Schneider et al.l I2010D has 
14 063 objects in this redshift range. 



G 



N. P. Ross et al. 



Description 



All DR9 boss_targetl quasar targets 6 
specPrimary = 1 

" and reliable redshift c 
" and 2.2 < z < 3.5 

XDQSO DR9 quasar targets 
" with spectra d 
" and reliable redshift c 
" and 2.2 < z < 3.5 
" and / BC > 0.85 



No. of Objects 
Pipeline" Visually 
Inspected 



182 973 
167 331 

132 290 163 128 
54 019 54 593 



74 607 
63 061 
54 416 
34 803 



62 048 
35 099 
23 301 



TABLE 2 

Properties of the SDSS-III BOSS DR9 QLF dataset. 
"The autom ated redshift determ ination algorithms 
described in iboltqn et al.i 120121) ■ '' the dr9 quasar 

boss.targetl TARGET FLAG IS DEFINED IN IRpSS ET AL.I ([2012). 
c TOTALS INCLUDE STARS. d ALL XDQSO DR9 QUASAR TARGET 
SPECTRA HAVE Specprimary=l BY DESIGN. THE UNIFORM SAMPLE 
DEFINED IN SECTION 2 IS BASED UPON THE XDQSO SELECTION. 



2.3. DR9 Uniform Sample 

We now define a uniform subsample from the parent 
DR9 quasar dataset. XDQSO models the distribution 
in SDSS flux space of stars and quasars as a function of 
redshift, as a sum of Gaussians convolved with photo- 
metric measurement errors, allowing the Bayesian prob- 
ability that any given object is a quasar to be calcu- 
lated. XDQSO is specifically trained and designed to 
select quasars in the redshift range 2.2 < z < 3.5 down 
to the BOSS limiting magnitude. 

XDQSO was only chosen as the algorithm to define 
the uniform sample after the first year of BOSS spectro- 
scopic observations. Each object is assigned a probabil- 
ity, P(QSOMIDZ), that it is a quasar with 2.2 <z< 3.5 
Objects with P(QSOMIDZ) > 0.424 ([Bow et al.ll20lll) 
are targeted as part of the uniform (CORE) sample. 
Knowing this threshold, we are able to say which targets 
XDQSO would have targeted in the first year of observa- 
tions, many of which BOSS did obtain spectra for. There 
are 74 607 quasar targets selected by XDQSO over the 
DR9 footprint, 63 061 of which have spectroscopic ob- 
servations, and of these, over half (35 099) ar e confirmed 
2.20 < z < 3.50 quasars by visual inspection ([Paris et al.l 

H3J). 

This sample of 35 099 quasars is over an order of 
magnitude more objects in this redshift range than 
in the study of iRichards et al.l (|2006bD from DR3 
(|Abazaiian et a l. 2005) . Fig. [2] shows the redshift distri- 
bution of this sample. Although we plot the full redshift 
range of the quasars, we do not use data from quasars 
which have a redshift below 2.2 or above 3.5, where the 
mid-z XDQSO selection, by design, is quite incomplete. 
The BOSS DR9 uniform quasar sample has a mean (me- 
dian) redshift of (z) = 2.59 (2.49). 

This uniform dataset is our primary basis for the QLF 
measurement. We supplement these data with a com- 
plementary dataset, selected by photometric variability 
criteria. 

2.4. Variability selection: Stripe 82 

The ~300 deg 2 area centered on the Celestial Equator 
in the Southern Galactic Cap, commonly referred to as 



"Stripe 82", was imaged repeatedly b y the SDSS over 
10 ye ars, generating up to 80 epochs (jAbazaiian et al.l 
2009), due in large part to the SDSS Supernova Sur- 
vey (jFrieman et al.1 12008T ) . These data are beneficial 
for quasar target selection for two reasons: (i) the im- 
proved photome try of the deeper data better defines 
the stellar locus (jlvezic et al.ll2007l ) an d fii) quasars can - 
be selected based on their variability (|Sesar et all 120071: 
i Bramich et alj|2008t [Schmidt et aLll2010HMacLeod et al l 
boilHPalanaue-Delabrouille et al.ll20lU l2012ri~ 

iPalanque-Delabrouille et al.l ([20111 ) describe the spec- 
troscopic quasar target selection for BOSS on 220 deg 2 of 
Stripe 82 based on variability. This variability selection 
was designed to select quasa rs with i > 18.0 and q < 22.3 
mag and redshift z > 2.15 (Palanau e-Delabrouille et al.l 
120111 Sec. 3.2). This dataset — which we shall re- 
fer to as the Stripe 82 (S82) data in what follows - 
includes ~6000 z > 2 quasars, roughly half of which 
would have been selected by XDQSO (as seen in Fig. [T]). 
Since the completeness of the variability selection is only 
very weakly dependent on redsh ift (e.g., Fig. 11 of 
IPalanque-Delabrouille et al"1l2011l) these data are subject 
to different, and arguably much weaker, selection biases 
than a color-based selection, as we show in Section [3T2l 

3. SURVEY COMPLETENESS 

To measure the QLF we must quantify the proba- 
bility, P(z, M, SED), of spectroscopically confirming a 
quasar of a given redshift z, absolute magnitude M 
and Spectral Energy Distribution (SED) shape. In this 
section, we describe the sources of incompleteness in 
the sample, and our checks of our completeness correc- 
tions. Our focus in this section will be on the DR9 
data; discussion of the completeness for the S8 2 sam - 
ple can be found in IPalanque-Delabrouille et all ([201 If ) 
and lRoss et"aH ([2012D . 

3.1. Incompleteness Descriptions 

We follow t he appr oaches of iCroom et al.l ([2004D and 
iCroom et al.l (|2009bD to quantify four avenues of po- 
tential sample incompleteness: Morphological, Target- 
ing, Coverage and Spectroscopic, and give descriptions 
of each type. 

3.1.1. Morphological completeness 

The input catalog to the BOSS quasar targeting al- 
gorithm is restricted to objects with stellar morpholo- 
gies in the single-epoch SDSS imaging. Host galaxies of 
z > 2 quasars are highly unlikely to be detected in this 
imaging, thus viable quasar targets should be unresolved. 
However, any true quasars not targeted because they are 
erroneously classified as resolved in the photometry will 
contribute to the survey incompleteness; this is referred 
to as morphological completeness (/ m ). 

We checked the assumption that z > 2 host galaxies 
are undetected in SDSS imaging, and tested the reliabil- 
ity of the star /galax y classifier from the SPSS photomet- 
ric p ipeline photo ([Lupton et al.l 120011 : iScranton et al.l 
12002( 1 . to the BOSS target selection magnitude limit. To 
do this, we compare to the Hubble Space Telescope (HST ) 
observations of the COSMOS field ([Scoville et alJl2007al) . 
which has been observed by the SDSS imaging camera. 
Objects classified as extended by SDSS that are actually 



The SDSS-III BOSS: Quasar Luminosity Function from DR9 



7 



TOT 

(J2000) 



~DicT 

(J2000) 



-band 



'- band 



r- band 



r band 



- band 



«- band 
extinction 



3 



^-pipc 



20.845±0.060 
22.591±0.295 
22.357±0.190 
21.434±0.139 
21.568±0.144 
21.097±0.108 
24.714±0.886 
19.307±0.028 
21.532±0.141 
22.735±0.319 



20.319±0.027 
20.455±0.029 
19.852±0.024 
20.876±0.030 
20.848±0.036 
20.439±0.027 
21.823±0.056 
18.788±0.020 
21.006±0.040 
21.830±0.057 



20.377±0.028 
20.013±0.026 
19.240±0.017 
20.805±0.038 
20.903±0.044 
20.349±0.032 
21.257±0.046 
18.736±0.018 
20.889±0.044 
21.918±0.082 



20.206±0.035 
19.686±0.030 
19.129±0.018 
20.648±0.042 
20.811±0.059 
20.216±0.038 
21.499±0.080 
18.747±0.022 
20.910±0.065 
21.769±0.096 



19.922±0.092 
19.650±0.081 
18.898±0.035 
20.214±0.132 
20.417±0.151 
19.772±0.090 
20.850±0.179 
18.571±0.034 
20.428±0.156 
21.639±0.275 



0.4615 
0.8947 
0.8333 
0.4615 
0.8947 
0.7143 
0.7143 
1.0000 
0.8947 
0.8947 



0.031620 
0.058656 
0.063211 
0.074886 
0.075538 
0.077683 
0.085803 
0.112584 
0.113820 
0.132704 



0.495352 
1.497665 
0.809249 
0.407500 
1.610326 
3.548377 
3.399193 
3.120975 
1.523919 
1.685750 



0.0527 
0.0509 
0.0585 
0.0540 
0.0485 
0.0563 
0.0569 
0.0451 
0.0500 
0.0476 



2.260 
3.228 
3.028 
2.282 
2.400 
2.237 
2.904 
2.353 
2.589 
2.526 



2.254 
3.228 
3.028 
2.281 
2.400 
2.238 
2.903 
2.343 
2.589 
2.526 



TABLE 3 

The BOSS DR9 statistical quasar dataset. The first ten lines are shown here for guidance regarding its format and 
content. The full table is published in the electronic edition of The Astro-physical Journ al. More details of the pipeline and 
visual inspection redshifts are documented in B olton et alJ Q2012I) and [Paris et al] 1120121) . 
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Fig. 3. — Splitting the sample of variability selected 2.2 < z < 3.5 Stripe 82 quasars into the 2 333 that are selected by the XDQSO 
algorithm (red) and those (3 143) that are not (blue). (Left): The distributions in the (u — g) vs. (g — r) color-color plane, (the stellar locus 
is given by the black contours); (Center): the (g — i) vs. i-band and (Right): the resulting N(z) histogram, (with all z_conf _person>3 
objects indicating a secure, visually inspected, redshift) plotted by the black line. 



unresolved in the COSMOS imaging, could be true high- 
z quasars that we fail to target in BOSS. We found that 
at r < 21.0, <3% of objects classified morphologically 
as galaxies by SDSS are unresolved in COSMOS; this 
fraction rises to f»8% at r = 22.0. We also found that 
all BOSS quasars at z > 2 lying within the COSMOS 
field are unresolved by HST. Thus we conclude that host 
galaxy contribution to morphological incompleteness is 
minimal, and we do not account for the misclassification 
rate of stellar objects by photo in our QLF calculations. 

3.1.2. Targeting completeness 

Targeting completeness, ft, accounts for any true 
quasars which are not targeted by our selection algo- 
rithm. We use the XDQSO method to select our high-z 
quasar targets, but, as we demonstrate below, the com- 
pleteness of XDQSO is a strong function of color, redshift 
and magnitude. Also, as we have noted, the XDQSO 
method was not used to select a uniform sample until 
the end of Year One. 

The area targeted with Year One target selection was 
1661 deg 2 , although the 220 deg 2 of Stripe 82 was re- 
targeted and re-observed in Year Two. iRoss et ail (|2012[ ) 
found that apart from over the Stripe 82 area, the frac- 
tion of objects selected by the XDQSO CORE algorithm 



which actually were targeted, was 87% for the DR9 foot- 
print. This result is consistent with the numbers of 
XDQSO targets (74,607) that have spectra (63,074), as 
given in Table [2] 

In Stripe 82, this fraction declines to 65.4%. This drop 
in targeting completeness is due to the deeper Stripe 82 
photometry which eliminates many noisy stellar contam- 
inants in the single-epoch XDQSQ target list, while se- 
lecting nearly all of the true quasars selected by CORE. 
The high targeting completeness fraction of XDQSO in 
the remainder of Year One is because many of the CORE 
quasar targets (and consequently true quasars) were se- 
lected by other target selection methods. There are 
in so me sense the "ea siest" quasars to discover. In- 
deed, IBovv et al.1 (|2011l) demonstrate tha t XDQSO and 
the Likelihood method (Kirkpat rick et aLll2011[ ). used as 
CORE for Year One, select similar samples. 

3.1.3. Coverage completeness 

Coverage completeness, f c (&), is defined as the fraction 
of BOSS quasar targets that have spectroscopic observa- 
tions, is quantified on a sector-by-sector basis, and is thus 
a function of angular position, 8. The main source of cov- 
erage incompleteness is fiber collisions, i.e. fibers cannot 
be placed closer than 62" to each other on a single plate. 
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Fig. 4. — The redshift distribution for radio-selected objects 
( "FIRST" , blue histogram) and those that passed both the radio 
and XDQSO selection ( "FIRST" +XD, green histogram). 



Description 


No. of Objects 




Pipeline 


V.I. 


All Stripe 82 quasar targets" 


15 576 




with specPrimary = 1 


12 576 




AND reliable redshifts 


10 506 


11 990 


AND 2.20 < z < 3.50 


5 433 


5 476 


AND w/ XDQSO scln. 


2 318 


2 333 



TABLE 4 

P roperties of the Stripe 82 BOSS QLF dataset, described 
in Palanque-Delabrouille et AL. (2011) AND WITH data from 
two of the Ancillary Programs that targeted quasars 
due to their near-infrared colors or radio properties. 



In Year One, the coverage completeness was > 90%, and 
in Year Two CORE quasars were given highest tiling pri- 
ority, and f c {9) is > 98%. 

3.1.4. Spectroscopic completeness 

Spectroscopic completeness, / s (mag, z, 9), is the frac- 
tion of BOSS quasar targets with spectra, from CORE, 
that have reliable redshifts. With the visual inspec- 
tions of all the quasar target spectra, this fraction is 
> 90%. We define a "spectro-coverage completeness" 
a s fsc = fc x fs, which is the fraction of XDQSO targets 
that were allocated fibers, and returned a reliable spec- 
trum. Tests showed that the computation of the QLF 
is only very weakly sensitive to the value of / sc , and we 
choose a threshold of / sc — 0.85 as a good compromise 
between high completeness and large sample size. Sec- 
tors with / sc > 85% are shown in red shades in Fig.[T] we 
limit our LF analysis to this area. This approach tends 
to exclude regions that have Year One spectroscopy, leav- 
ing an area of 2236 deg 2 . There are 23 301 quasars 
in this area (all with visually confirmed redshifts) 
and this sample is given in Table [3j 

3.2. Empirical checks using Variability selected 
Quasars 

The 220 deg 2 of spectroscopy across the Stripe 82 
field has targets selected via the ir optical variability 
(|Palanque-Delabrouille et al.|[201lT and g[2Ti|). We con- 
centrate on the 5 476 2.20 < z < 3.50 quasars in 
Stripe 82 (Tabled]), including 122 quasars selected solely 
due to their near-i nfrared colors or radio properties 
(jDawson et al.ll2012D . In this area, we find a higher sur- 
face density of high-z quasars, 24.9 deg" 2 , than across 
the full DR9 dataset (14.7 deg" 2 ) and the XDQSO uni- 
form sample (9.6 deg -2 ). Thus, this enhanced Stripe 



82 dataset is more complete and less affected by color- 
induced selection biases, and we will use it to measure 
the targeting completeness of our XDQSO uniform sam- 
ple empirically. 

We split the sample of 5 476 visually confirmed 2.2 < 
z < 3.5 quasars into the 2 333 that would have been 
selected by the XDQSO algorithm (XD) and those (3 143) 
that would not have been (!XD). Over 96% of the !XD 
sample was selected by a variability algorithm. The (u — 
g) vs. (g—r) color-color plane, their distribution in (g—i) 
vs. i-band and the resulting N(z) redshift histograms of 
these two samples are given in Fig. [3] 

The difference between t he two selections is apparent 
and co nsistent with that in Palanqu e-Delabrouille et all 
(|2011l their Fig. 18). The !XD sample more heavily 
overlaps with the stellar locus in (u — g) vs. (g — r), 
and is generally redder than the XD population in (g — r) . 
Thus, the variability selection is able to recover quasars 
from the stellar locus. The distribution of XD and !XD 
objects in the (g — r) vs. (r — i) color-color plane (not 
shown) is similar, in that the ! XD population overlaps 
with the stellar locus, but both populations have similar 
distributio ns in the (r — i) color (again in agre ement with 
Fig. 18 of Palanqu e-Delabrouille et afl 120111 ). However, 
the !XD population is redder in (g — i) (Fig. [31 center). 

The N(z) histograms for the two sampl es are also very 
much i n line with previous studies (e.g., Richards et al.l 
I2006bt iPalanque-Delabrouille et al.l l201lf T The decre- 
ment at z ~ 2.7 — 2.9 for the XD selection is due to the 
overlap of such objects with the stellar locus in color- 
space, a key issue in the original studies in SDSS. 

The two samples are similar in their distributions 
across z = 3.0 — 3.3, though at z m 3.4 — 3.5, there 
are more !XD objects, probably due to the efficient cut- 
off of the z = 1.2- 3.5 "mid-z" XDQSO selection, and 
the fact that at z ~ 3.5 quasar colors again approach the 
stellar locus. 

3.3. Radio Selection versus Color Selection 

Figure [4] displays the redshift distribution for radio- 
selected objects ("FIRST", blue histogram) and those 
that passed both the radio and XDQSO selection 
( "FIRST" +XD, green histogram). This graph can be 
compared directly to Figure 10 in R06, which com- 
pares the redshift distribution of radio-selected quasars 
to those that were both radio- and colo r-selected us- 
ing the full DR3Q ([Schneider et al.l 120051 ) sample. The 
redshift distribution of the radio-only selected objects 
is smoother and has a smaller decrement of objects at 
z = 2.7 — 2.8 than the radio+color selection. This was 
also seen in the R06 DR3Q investigation. 

However, we are wary of over-interpreting this for sev- 
eral reasons Firs t, only ~ 2% (3 348) of the DR9Q 
([Paris et al.l HoTl), and - 2% (747) of the XDQSO 
quasars are targeted via their radio properties, of which 
half are selected only via their radio properties. Second, 
BOSS is deeper than SDSS, whereas the FIRST detection 
limits are the same for the two optical surveys, so BOSS 
radio sources are more radio loud. If radio l oudness cor- 
relates with redshift and/or lu minosity (e.g.. Uiang et al.l 
l2007tlSingal et al.ll201lLl2012l ). an attempt to correct the 
N(z) distribution using radio-loud quasars would be in- 
correct (see also the discussion in § 3.4 of R06). Finally, 
radio-loud quasars are not drawn from the same color 
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Fig. 5. — Color-redshift relation for BOSS quasars selected by XDQSO (solid, blue) and for simulated quasars (dashed, red). The thick 
lines show the median values, and the thin lines show the 25%- 75% range, at each redshift. The agreement between the data and models 
is generally very impressive. At z < 2.2 the colors are less well matched, but are very noisy due to the very low completeness of XDQSO 
at those redshifts. 



distribution as radio-quiet quasars, with t he radio-loud 
popu l ation tending to have redder colors (iWhite et al.1 
l200l IMcGreer et alJ 120091 : IKimball et afl 120 1 ID . Thus 
objects in a radio-selected quasar sample do not share the 
same selection function as a purely color-selected sample. 

3.4. Simulated Quasar Spectra and Completeness 

In Section [3~2l we presented the sample established by 
the XDQSO targeting algorithm, and compared it to a 
dataset constructed from the Stripe 82 sample of quasars 
selected independently of that algorithm. We can use the 
results of § 13.21 to quantify the completeness of XDQSO 
only in the limit that the Stripe 82 sample is itself com- 
plete. Here we adopt another approach: we construct 
a model for the observed spectroscopic and photometric 
properties of quasars, generate a large sample of sim- 
ulated quasars, and then test the targeting algorithm 
against this model using the simulated quasars (e.g-. lFanl 
S IRichards et alJl2006bh . 

The broadband optical fluxes used by XDQSO are 
dominated by a featureless power-law continuum. How- 
ever, quasar selection is highly sensitive to the colors of 
quasars, which evolve strongly with redshift as the broad, 
high-equivalent- widt h emission lines move through the 
optical bandpasses (|Richards et al.l l2001a[ ). Hence, a 
complete prescription for quasar properties must capture 
both the smooth continuum and the emission lines. 

Past models have generally adopted a continuum power 
law index of a v = —0.5, where a v is the frequency 
power law index (i.e. F(u) oc y 01 "), typically mea- 
sured from quasar spectra fe.g.. |Richstone fc Schmidt! 
[19801: IFrancis et al.lll99lt IVanden Berk et al.H200lD . Of- 
ten a break is added to the near UV where a softer 
spectrum is observed (a ny ~ -1.7, iTelfer et all [20021 : 
iShang et al.ll2005l 1201 lh . Emission line templates, in- 
cluding Fen complexes, are then generated from com- 
posite me an quasar spec tra c onstructed from larg e sam- 
ples (e.g., IFrancis et aT1ll99lH\^nden Berk et alj|200lD . 
To these emission components, absorption from the Lya 



forest is added, given a model for its redshift dependence. 
With these basic assumptions, models can be generated 
that broadly reprod uce the mean colors o f quasars as a 
function of redshift (|Richards et al.ll2001al ). 

For this work we have taken advantage of the 
many improvements in our understanding of quasar 
spectral properties in recent years, namely, improved 
meas urements of absorption due to the Lya forest 
(e.g., iWorseck fc P rochaska 120 111 ), templates for iron 
emission (a significant co ntributor to quasar colors; 
iVestergaard fc Wilkes 2001) and finally, large samples of 
quasars to calibrate the models. We have simulated the 
full survey, by passing the model quasars through the 
target selection algorithm and comparing the resulting 
color distribution to observations. Under the common 
assumption that quasar spectral features do not evolve 
with redshift, the selection function provides a redshift- 
dependent window into the underlying color distribution. 
By comparing the colors of quasars that the model pre- 
dicts are selected by the survey, to those actually ob- 
served, we can determine a best-fit model that not only 
recovers the selection function, but also provides insight 
into the intrinsic properties of quasars. 

This process will be detailed in a forthcoming work 
(McGreer et al., in prep). Here we briefly outline the 
steps taken to generate a model for the population of 
quasars observed (and not observed) by BOSS. 

1. We construct a grid of model quasars in 
(M ,z) space, u sing t he luminosity function from 
iHopkins et all ([20071 ) . For each quasar we ran- 
domly sample the following components: 

(a) A broken power-law continuum with a break 
at HOOA; at near-UV wavelengths the power 
law index is drawn from a Gaussian distri- 
bution with mean a u = —1.7 and scatter 
a{a v ) = 0.3; for A > 1100 A the distri- 
bution has a mean a v = —0.5 and scatter 
cr(a„) = 0.25. 
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Fig. 6.— The selection function for the BOSS XDQSO sample 
via simulated quasar spectra and photometry. Contour levels are 
drawn at 1, 5, 10, 20, 40, 50, 60, 80, and 90 percent complete- 
ness, as determined by the fraction of simulated quasars selected 
by XDQSO as a function of redshift and i-band magnitude. The 
50% and 90% levels are drawn with dashed lines. 



(b) Emission lines generated from composite spec- 
tra of BOSS quasars binned in luminos- 
ity, reproducing trends between emission line 
properties and con tinuum luminosity (e.g., 
the Baldwin Effect: iBaldwinl fl977t IWu et alJ 
200§). The resulting emission line template 
provides the mean and scatter in line strength 
for prominent quasar emission lines as a func- 
tion of luminosity; values for individual ob- 
jects are drawn from this distribution. 

(c) Fe emission from the template of 
iVestergaard fe Wilkesl (|200l . The tem- 
plate is divide d into discrete wavelength 
segments (see IVestergaard fe Wilkesl [2001) 
that are scaled independently; the scale 
values are determined during the fitting of 
the composite spectra used for the emission 
line template. 

(d) Lya forest blanketing according to the pre- 
scription of iWorseck fe Prochaskal (|2011[ ). A 
population of absorbers is generated in a 
Monte Carlo fashion using th e para meters 
given m IWorseck fe Prochaskal (|2011f l. The 
lines are modeled as Voigt profil e s usin g the 
approximation of Te pper-Garrial ()2006f ). and 
then applied to the forest regions of the sim- 
ulated spectra. All Lyman series transitions 
up to n = 32 are included. A total of 5000 in- 
dependent sightlines were generated and then 
randomly drawn for each of the simulated 
quasars. 

2. We generate spectra from this grid and calculate 
SDSS broadband fluxes from the spectra. 

3. The fluxes are transferred to observed values via 



empirical relations for the photometric uncertain- 
ties derived from single-epoch observat ions of stars 
on St ripe 82, using the coadded fluxes ([Annis et al.1 
120111 ) as the reference system. 

4. The XDQSO algorithm is used to calculate mid-z 
quasar probabilities for each model quasar in the 
same manner as for BOSS selection, and a sample 
of model quasars is defined. 

This describes our fiducial model. We further test two 
modifications to the fiducial model. For comparison to 
previous work, we implement a second model where the 
emission line template is derived from a single compos- 
ite spectrum and thus does not have any dependence on 
luminosity. This template c omes from the SPS S com - 
posite spectrum presented in iVanden Berk et al.1 (j2001[ ) 
and is referre d to as "VdB l i nes". This model is clos- 
est to that of Ric hards et al.l (|200 6hV Finally, we test a 
third model that includes dust extinction from the host 
galaxy. In this model, in dividual spectra ar e extincted 
using a SMC dust model (jPrevot et al. 1984), with val- 
ues of E(B — V) distributed expo nentially around a peak 
of 0.03 (e.g.. lHopkins et al.l l2004'). This model is referred 
to as "exp dust" . We compare the three models in more 
detail in Appendix [A"l 

We test the accuracy of the fiducial model by check- 
ing the simulated quasar colors against observed quasar 
colors. This prior is used to distribute the simulated 
quasars in flux and redshift space in a manner similar to 
the intrinsic distribution. The simulated quasar photom- 
etry is passed through the XDQSO selection algorithm 
to mimic the observations, so that the final color rela- 
tions for simulated quasars are derived only for objects 
that would have been targeted by the s urvey. We then 
constr uct the color- redshift relation (e.g.. lRichards et all 
l2001a|) of both simulated and observed quasars by di- 
viding the samples into narrow redshift bins (Az = 0.05) 
and calculating both the median and scatter of the u — g, 
g — r, r — i, and i — z colors within each redshift bin. 
The results for the fiducial model are shown in Figure El 
demonstrating that the model does an excellent job of 
reproducing the observed quasar color distribution. 

Dust extinction is thought to produce the red tail of the 
color d istribution often seen in qu asar surveys. For ex - 
ample, iRichards et"aT1 (|2003| ) and iHopkins et all fj2004D 
find that ~ 20% of SDSS quasars have colors consistent 
with reddening from dust with an SMC-like extinction 
curve with E(B — V) > 0.1. We find that the exp dust 
model does not significantly improve the fit to the color 
distribution of BOSS quasars compared to our fiducial 
model, and thus our primary analysis does not include 
the effect of dust extinction. Section [5] will explore how 
the differing assumptions of the three models affect the 
calculation of the QLF. 

Tableland Fig. [5] give the selection function, i.e. the 
fraction of selected quasars, in each bin of Al and z, 
generated from the fiducial model outlined above. 

We compare the model selection function to an em- 
pirical relation from Stripe 82 in Figure The green 
lines show the fraction of model quasars that are selected 
by XDQSO-CORE. This is compared to the fraction of 
Stripe 82 quasars — predominantly selected by variabil- 
ity criteria — that are recovered by XDQSO selection. 
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Fig. 7. — BOSS (XDQSO) quasar selection function, in discrete redshift bins covering the range 2.2 < z < 3.4. The blue points with error 
bars show the empirical selection function derived from Stripe 82; specifically, they denote the fraction of Stripe 82 quasars selected by 
XDQSO within bins of magnitude and redshift (the x-error bar rep resen ts the i- magnitude bin width). The green line shows our "fiducial" 
model selection function using simulated quasars as described in § 13.41 , and the same binning as the empirical points. The agreement at 
z < 3 shows that the two are consistent, as expected if both the model is a good representation of true quasars and the variability selection is 
highly complete. The model predicts lower efficiency at z > 3, suggesting that the completeness of the variability selection is lower at higher 
redshifts i.e., XDQSO recovers a higher fraction of the variability quasars than the model quasars, and variability is potentially missing a 
population of quasars. Note also that the efficiency predicted by the model is generally lower in the faintest magnitu de bin (i ~ 22), again 
suggesting that variability was less complete at the faint end. Note that the model has a luminosity function prior (Hopkins ct al. 2007) 
applied. For comparison, we also plot the "VdB lines" and "exp dust" model selection functions (red and cyan lines, respectively), both of 
which generally show poorer agreement with the empirical points from Stripe 82 than the fiducial model. 
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TABLE 5 

The Quasar Selection Function for the fiducial model 
described in the text; see also flg. [6] the final column 
gives the fraction of simulated quasars selected by 
XDQSO. Table[5]is published in its entirety in the 
ELECTRONIC EDITION OF The Astrophysical Journal; THIS 
EXCERPT HERE IS SHOWN HERE FOR GUIDANCE REGARDING ITS 
FORMAT AND CONTENT. 



This empirical relation is shown as blue squares with er- 
ror bars (Poisson uncertainties) . The agreement at z < 3 
shows that the two are consistent, as expected if both the 
model is a good representation of actual quasars, and the 
variability selection is highly complete. 

However, there is some disagreement between the two 
completeness estimates. For example, smaller fractions 
of model quasars are selected by XDQSO at z > 3 than 
are selected by XDQSO from the Stripe 82 sample in the 
same redshift and magnitude bins. This may indicate a 
deficiency of the models; however, we are encouraged by 
the excellent agreement between the colors predicted by 
the model and those observed (Figure [5]). Alternatively, 



our assumption that the variability-selected sample is 
both complete and unbiased may be invalid. Indeed, 
color criteria were applied to objects when construct- 
ing t he variability sample ([Palanaue-Delabrouille et al.l 
1201 If ): these criteria may exclude some populations of 
quasars, in particular, they may introduce bias against 
high redshift quasars. In that case, the disagreement in 
Figure[7|suggests that XDQSO recovers a higher fraction 
of the variability-selected quasars, but both XDQSO and 
variability are missing a population of objects. This ef- 
fect may also explain why the model predicts lower com- 
pleteness at the faintest magnitudes: i.e., both XDQSO 
and variability have low completeness at i ~ 21.8, but 
XDQSO recovers a higher fraction of the quasars that 
are also selected by variability. 

In what follows, we implement our fiducial selection 
function model to calculate the QLF from the DR9 uni- 
form quasar sample. Since the color selection incomplete- 
ness dominates over the other sources of incompleteness, 
we do not make any further corrections during the QLF 
calculation. 

4. NUMBER COUNTS AND ^-CORRECTIONS 

4.1. Number Counts 

In Fig. [5] we present the cumulative i-band number 
counts of the datasets described in Sec. 2: the XDQSO 
uniform sample of 23 301 quasars across 2236 deg 2 with 
2.2 < z < 3.5 (red circles) and the 5 470 quasars 
across 220 deg 2 of Stripe 82 also with 2.2 < z < 3.5 
(red crosses). Also shown are the number counts from 
1.0 < z < 2.2 quasars selected from deeper, g » 23 spec- 
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Fig. 8. — Cumulative i-band number counts. Here, the 2.2 < 
z < 3.5 BOSS samples are in red, with the uncorrected BOSS 
uniform sample shown by the open circles, while the Stripe 82 
data are given by the crosses. Also shown are the number counts 
from the deeper, g ?s 23 1.0 < z < 2.2 quasars selected fro m the 
"boss21+MMT" survey UPalanqu c-Dclabr ouTlle et al.ll20"T3 ). For 
clarity, wc only plot errorbars at the bright end (i < 19), since the 
errors are smaller than the points at the faint end. We also show 
the double power law fits to the data as described in the text. 
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Fig. 9. — Comparison of the fc-correction from our fiducial model 
(blue line) with the R06 fe-correction (green line). Both are defined 
as the correction required to transfer the observed i-band flux to the 
i-band luminosity at z = 2, i.e., M; (z = 2). The offset of ~0.1 mag 
at z < 2.7 is due to a different treatment of Fe emission, and grows 
somewhat larger at higher redshifts as the C III] line enters the i- 
band. The blue line shows our quasar model for Mi(z = 2) = —26, 
and the shaded regions show the variation of the fc-correction with 
luminosity over the range —27 < Mi(z = 2) < —24.3, covering 
most of the luminosity range of BOSS quasars. Finally, the dashed 
line shows a pure continuum fc-correction for a v = —0.5. 



troscopy, using da ta from the "boss21+MMT" surv ey 
(blue, filled circles: IPalanque-Delabrouille et al.ll2012t ). 

For the uniform BOSS sample, the open red circles 
are for the raw, uncorrected number counts, whereas the 
filled circles use the correction derived in the previous 
section, integrated over our redshift range. These num- 
ber counts can be compared to the 2.2 < z < 3.5 quasars 



z cm fc-corrcction 

0.105 0T323 

0.115 0.250 

0.125 0.317 

0.135 0.332 

0.145 0.335 

0.155 0.334 

0.165 0.285 

0.175 0.291 

0.185 0.334 

0.195 0.249 



TABLE 6 

THE i-BAND fc-CORRECTIONS. THE fc-CORRECTION IS OBTAI NED 
USING THE FIDUCIAL QUASAR MODEL DESCRIBED IN SECTION 13.41 

AND INCLUDES AN UPDA TED TREATMENT O F THE EM ISSION LINE 
TEMPLATE COMPARED TO IRlCHARDS ET AlI ||2006bI) . We DEFINE 
OUR fc-CORRECTION TO BE OUR MODEL fc-CORRECTION AT 

Mi(z = 2) = —26.0 (see main text for details). The full 

TABLE IS PUBLISHED IN THE ELECTRONIC EDITION OF The 
Astrophysical Journal; THE FIRST TEN LINES ARE SHOWN HERE 
FOR GUIDANCE REGARDING ITS FORMAT AND CONTENT. 



selected via their variability signature on Stripe 82. The 
two are in reasonable agreement to i ~ 21.0, with the cor- 
rected number counts being consistently higher. Fainter 
than this, the variability number counts drop more no- 
ticeably below the corrected counts, suggesting that this 
dataset is incomplete at the faint end (or that the incom- 
pleteness of the DR9 sample is overestimated). Across 
the redshift range 2.2 < z < 3.5 and down to i = 21.5, 
the corrected BOSS DR9 uniform cumulative number 
counts reach 34.4 deg -2 , whereas the Stripe 82 cumu- 
lative counts are 26.2 deg -2 . 

Motivated by the double power-law form of the QLF 
(Eqn . [2]), and prior measurements (e.g. I M vers et al.l 
2003). we also express the cumulative number counts as 
a double power-law, 



J 10- Q <i( m - m o) -)- \§-Pd{m-m ) v ' 

and find best-fits to the (corrected) BOSS uniform and 
boss21+MMT counts. For the BOSS sample we find 
slopes of ad — 1-50 and fid = 0.40, while the "break 
magnitude" too =19.0 and the normalization, Nq = 2.63 
deg -2 . In comparison, the boss21+MMT data has a 
less-steep bright end slope of ad — 0.80 and an al- 
most flat faint end slope (3d = 0.10. The break mag- 
nitude is fainter at too = 20.4 and the normaliza- 
tion is significantly higher, Nq — 43.6 deg -2 . These 
power-law descriptions and surface densities will allow 
extrapolatio n for future Lyq-forest c osmology experi- 
ments (e.g., IMcQuinn fc White] 1201 ID . For unobscured 
1.0 < z < 2.2 quasars, there are 48 (78) objects deg -2 
down to i < 21.5 (23.0), broadly consistent with the 
value of 99±4 quasars deg -2 with gdered < 22.5 from 
iPalanqueT Mabro uille et al.l (|2012t their Table 5) and a 
surface density si milar to that selected by a shallow m id- 
infrared selection ()Stern et al.M2012t lYan et al.ll2012t) . 

4.2. k-corrections 

Following R06, we define the fc-correction to deter- 
mine the i-band luminosity at z — 2; Mi(z — 2). This 
is ^2700 A in the rest-frame, and is close to the me- 
dian redshift of the BOSS quasars sample. R06 calculate 
the quasar fc-correction as the sum of a component due 
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Fig. 10. — The Absolute Magnitude-Redshift distribution for 
SDSS DR7 quasars (black points), BOSS DR9 (red points) and 
the fainter, boss21+MMT variability-selected datasct from Stripe 
82 (blue squares). The normalized redshift distributions of the 
three datasets are shown in the bottom panel, while the normal- 
ized absolute i-band magnitudes are given in the side panel. The 
bright and faint magnitude limits of the BOSS DR9 sample are 
given by the solid turquoise lines; the wiggles are due to the red- 
shift dependence of the fc-correction. 



to the underlying continuum, k cont , and a component 
due to the emission lines, fc em - The sign convention of 
the fc-correction, k(z\ is rointnnsip = WnhssrvBri — k(z) 
(|Oke fe SandageHl968t iHogg et all 
We obtain a fc-correction from 



12001 . 

the fiducial 



quasar 

model defined in Section 13.41 This model is very sim- 
ilar to the one adopted by R06. The continuum model is 
identical: a power-law slope of a v = —0.5 at A > llOOA, 
and we set k cont (z = 2.0) = 0.0 by definition. On the 
other hand, our emission line template is not the same 
as the one used by R06. They defined their emission line 
template using a single composite sp ectrum derived from 
SDSS quasars (similar to that of Vanden Ber k et aLl 
2001) with a power-law continuum removed. Our emis- 
sion line template is similarly obtained from composite 
spectra; however, we use a suite of composite spectra 
binned in luminosi ty, and fit the continuum join tly with 
the Fe template of iVestergaard fe Wilkesl (|200lD . 

One advantage of defining the fc-correction to be in 
the i-band at z = 2 is that the i-band is relatively free of 
strong emission lines at this redshift. At z < 2.5, the i- 
band samples rest-frame ~ 2200A, where the only strong 
emission line features are from Fe II and Fe III. We find 
that our model for Fe emission introduces a shift of about 
~0.1 mag at 2 < z < 2.7 relative to the R06 model; i.e., 
fcfid — fcR06 ~ —0.1. At higher redshifts, the Cm] line 
enters the i-band and the offset between our fc-correction 
and that of R06 grows somewhat larger, reaching « — 0.2 
mag at z = 3.5. We compare our fc-correction to that of 
R06 in Figure [9] 

Though our quasar model introduces a luminosity de- 



pendence to the fc-correction due to the anticorrelation 
between emission line equivalent width and luminosity 
(the Baldwin Effect), we chose not to apply this fur- 
ther correction in this work. For the same reasons given 
above, at z < 2.7 there is almost no variation in our fc- 
correction with luminosity. At z = 3.5, this effect only 
reaches ~ 5% over the range of luminosities probed by 
BOSS (see Figure[9]). This variation is much smaller than 
the intrinsic scatter in fc-corrections at a given redshift; 
i.e., even if we attempted to correct for the Baldwin Ef- 
fect in the mean, the scatter in this correlation is far 
greater than the correction. We define our fc-correction 
to be our model fc-correction at Mi(z = 2) = —26.0, 
near the median luminosity of the BOSS sample. Ta- 
ble [6] presents our new fc-correction. 

5. LUMINOSITY FUNCTIONS 

In Fig. [TUl we show the coverage in the ab- 
solute magnitude-redshift (M< — z) plane for the 
three datasets of main interest h er e: the SPSS 
(blac k points; I Richards et al.l l2006bl : ISchneider et al.l 
120101) : the XDQSO-selected 2.2 < z < 3.5 BOSS 
DR9 sample (red points) and the fainter variability- 
selected d ataset from the boss21+MMT sa mple (blue 
squares; IPalanque-Delabrouille et al.l 120121 ). We also 
analyze the Stripe 82 var i ability -selected dataset of 
IPalanque-Delabrouille et al.l (|201ll ). which has a similar 
redshift distribution as the DR9 sample. The bright and 
faint magnitude limits of the BOSS DR9 sample, i = 17.8 
and i — 21.8, respectively, are given are given by the solid 
turquoise lines. Our binning is identical to R06; the edges 
of the redshift bins in which we will calculated the QLF 
are: 0.30, 0.68, 1.06, 1.44, 1.82, 2.20, 2.6, 3.0, 3.5, 4.0, 
4.5, and 5.0, and the Mi bins start at -22.5 and are in 
increments of 0.30 mag 28 . 

5.1. The Optical Luminosity Function to i — 21.8 

In Fig. [TT]we show the i-band luminosity function from 
our BOSS DR9 uniform sample over 2.2 < z < 3.5, as 
well as the fainter boss21+ MMT sample of quasars cov- 
ering 0.68 < z < 2.2 from Palanqu e-Delabrouille et all 
(120121). We use th e binned QLF estimator 29 of 
IPaee fe Carreral (l2000l) . 



N„ 



(dV/dz) dz dL 



(4) 



This involves calculating the number of quasars, N q ob- 
served in a given (Mi — z) bin, correcting for our selection 
function, and dividing N q by the effective volume element 
dV of that bin. The effective volume is calculated by 
using our fixed, flat (fi A ,fi m ,/i) = (0.70,0.30,0.70) cos- 
mology, and the area of our uniform DR9 sample (2236 
deg 2 ). We check, and find that our redshift bins are suf- 
ficiently narrow to avoid complications due to evolution. 
The plotted error is estimated by 



SN 



"est 



dz dL 

28 All the necessary data and code used here to produce our 
results will be publicly available at www . sdss3 . org/dr9/qlf 

29 Croom et al. (2009a) show that the difference between the 
Page & Carrcra (2000) estimator and the "model- weighted" esti- 
mator of IMivaii et al.l (120011 ) is small, even at the bright end, for 
the z > 1 QLF. 



(5) 
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Fig. 11. — The i-band Quasar Luminosit y Function. The red points are from our analysis of B OSS quasars from DR9, while th e black 
squares are from the DR3 analysis of Richards ct al. (2006b). The boss21+MMT sample from Palanque-Dclabrouille et al. (2012) is also 
shown for 0.68 < z < 2.20 (blue filled squares). Over the redshift range 2.20 < z < 3.50, we use the 23 301 DR9 quasars uniformly selected 
by XDQSO, and that are in sectors of spectro-completeness of 85% or higher. The solid line in each panel is the BOSS DR9 QLF at 
2.2 < z < 2.6, to show how the luminosity function evolves. The open circles show the 2.2 < z < 3.5 QLF without correcting for the 
(fiducial) selection function. There are no uniform DR9 measurements above z = 3.5, since the XDQSO selection deliberately cuts off at 
this redshift. The Poisson error bars for the BOSS measurements in the three panels spanning 2.2 < z < 3.5 are the same size, or smaller, 
than the points shown. 



and SN is given by Poisson statistics including the up- 
weighting by the inverse of the completeness. We discuss 
the validity of this error estimate below. The binned 
QLF is also given in Table which gives the mean red- 
shift of the quasars in each bin, the mean i-band mag- 
nitude of the quasars in the bin, the magnitude at the 
bin center, the raw number of quasars in the bin, the 
log of the space density, $, in Mpc~ 3 mag -1 , and the 
error xlO 9 . The results from R06 using t he SPSS DR3 
are gi ven as the black squares in Fig. [TT] IShen fc Kellvl 
pOia) measured the QLF from the final DR7 SDSS 
quasar sample, and found excellent agreement with the 
DR3 results. 

Where the surveys overlap, we generally see very good 
agreement with the BOSS and SDSS data points, espe- 
cially at z < 3. Although there is overlap in L — z cov- 
erage between the SDSS DR3 and BOSS measurements, 
since DR3 and DR9 cover different ares of the sky, there 



are only 304 quasars (< 2%) are common to both sur- 
veys, mostly in the 3 < z < 3.5 redshift range. 

The limiting magnitude for BOSS DR9 quasar targets 
is g < 22.00 or r < 21.85, and with (r - i) « 0.05 for 
quasars at z ~ 2.5, we show an z-band limiting of i = 21.8 
as a guide in Fig. [TTJ There is strong evidence for a turn- 
over in the QLF, well before this limit, seen in all the red- 
shift panels i.e. up to z — 3.5. We shall see in Sec. 15.21 
that our results across 3.0 < z < 3.5 are also consis- 



tent with a turn over seen in other experiments, 
turn-over has been seen in the X-ravs (Mivaii et al.l 


This 
2001; 


Ueda et al. 


2003 


, Hasingcr ct al. 2005: Aird et al. 


2010; 


Fiore et al. 


2012 


) and also in the optical (Bovle et al. 



19881; ICroonTe t al. 2004, 2009a), with the BOSS DR9 



now extending this evidence to redshifts z = 3.5. 

Our calculation of the errorbars given in Fig. [11] as- 
sumes the errors in each bin are independent and are 
dominated by Poisson statistics of the observed objects. 
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Fig. 12. — The sensitivity of the BOSS DR9 Quasar Luminosity Function to the selection function models used. The QLF is divided by 
the fiducial model. The red points show this fiducial model (at by definition), with the extent of the points representing the statistical 
uncertainty range. The QLFs derived from the other two models are divided by the fiducial model QLF to highlight the effect of the choice 
of selection function on the derived QLF. The VdB lines model is shown as green points, and the exp dust model as blue points. In general 
there is good agreement between the three models, but in some redshift bins the disagreement can be 20% or greater. The fiducial model 
provides the best fit to the observed color-rcdshift relation; however, the differences seen here quantify the systematic uncertainty inherent 
in not knowing the selection function exactly. 



For the DR9 sample overall this is reasonable, given the 
very large volume surveyed (which reduces fluctuations 
due to large-scale structure) and the low mean occupancy 
of quasars in halos (which reduces the impact of halo 
count fluctuations on the correlations). When comparing 
to surveys of smaller volume, sample variance may dom- 
inate over the Poisson errors. In some redshift ranges, 
however, BOSS is quite incomplete, and require a signifi- 
cant selection function correction (compare the open red 
circles to the filled red circles in the 2.6 < z < 3.0 bin of 
Fig. [TT] for example). In these bins the error is dominated 
not by Poisson statistics but by the uncertainty in our 
estimate of the selection function (see Fig. [7]). This un- 
certainty can reach 50%, fractionally, for faint quasars in 
the most incomplete redshift range, leading to a similar 
fractional uncertainty in the QLF. However for most of 
the range plotted the uncertainty is significantly smaller. 

In Fig. Q/2] the effect of the selection function correc- 
tion is investigated further. Here we plot the logarithm 



of the ratio of the QLF number densities for the two other 
selection function models, "VdB lines" and "exp dust", 
introduced in Sec. 13.41 compared to our fiducial model 
(that is used to calculate the QLF presented in Fig. [TTJ. 
We concentrate on the redshift range 2.20 < z < 3.50. 
The errorbars for each model represent the Poisson un- 
certainties, and the differences between selection function 
models dominate over these statistical uncertainties, es- 
pecially at the faint end. The corrections derived from 
the exp dust model generally augment the estimated lu- 
minosity function, particularly at low luminosities and 
higher redshifts. This is likely due to the fact that BOSS 
quasar selection is flux-limited in the g and r bands, so 
that fainter and higher redshift objects subjected to dust 
reddening will be extincted out of survey selection. The 
corrections derived from the VdB lines model show an 
even stronger trend with luminosity. The dependence of 
observed quasar colors on intrinsic luminosity resulting 
from the Baldwin Effect leads to a luminosity-dependent 
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Fig. 13. — The i-band Quasar Luminosity Function. Black squares are from Richards ct al. (2006b), the rilled red circles are our 
measurements from the full DR9 sample, while the teal points are from our analysis of t he 5 ,476 BOSS quasars with 2.2 < z < 3.5 located 
and observed on Stripe 82. The top-left panel is the 2.2< z <2.6 measurement from Fig. [TT] shown as a guide. 



selection function. A QLF estimate that does not ac- 
count for this effect will incur an artificial tilt as a func- 
tion of luminosity, as highlighted by the figure. This tilt 
will further affect QLF parameters such as the power law 
slopes. 

In Fig. [T31 we continue to concentrate on the redshift 
range 2.2 < z < 3.5, and divide it more finely in redshift 
than in Fig. QTJ The data displayed in Fig. [13] are pre- 
sented in tabular form in Appendix[Bj and it will be these 
data that we will fit models to in Sec. 16.11 We compare 
the BOSS DR9 measurements to the 5 476 2.2 < z < 3.5 
quasars on Stripe 82 that were selected via variability. 
The BOSS DR9 and Stripe 82 measurements are in very 
good agreement below z ~ 2.7, consistent with the selec- 
tion function agreement in Fig. [7] 

However, there are differences between the two 
datasets for z > 2.7, especially at the fainter end. The 
DR9 measurement implies a higher space density than 
the Stripe 82 variability measurements. One possible 
explanation could be that the selection function is un- 
derestimated (in the sense that it over corrects N q ) from 
Sec. 13.41 However, this would potentially lead to higher 
DR9 space densities at the faint end at all redshifts. An- 



other possibility is that the variability selection is begin- 
ning to break down at the faint end, as the selection is 
based on light-curves taken from single-epoch imaging, 
and is susceptible to imaging incompleteness. 



5.2. Comparison to Other Results 

In Fig. [H we compare our BOSS DR9 QLF to other 
measurements of the QLF at z > 2. In each panel, we 
divide the QLFs by our best-fit "log-linear" Luminos- 
ity Evolution and Density Evolution (LEDE) model, de- 
scribed in Section 16.11 We concentrate on the redshifts 
z ~ 2. and z « 2.4, and the results of ICroom et al.l 
(|2009aD from the 2SLAQ QSO survey. We also compare 
our measurements a t z s=s 3.2 with recent results from 
IMasters et ail (12012ft using observations from COSMOS 
HScoville et aLll2007bft . We additionally compare the re- 
sults f rom our sister study, IPalanque-Delabrouille et al.l 
(|2012ft . at all three epochs. 

The 2SLAQ results are presented in ICroom et al.l 
(|2009aft a s a function of MJz , = 2), and the COSMOS 
results in IMasters et al.l (|2012ft in M1450, so in order to 
make a direct comparison, we convert these to Mi(z = 2), 
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Fig. 14. — The BOSS DR9 Quasar Luminosity Function compare d to other surveys. L eft: Measurement in the 1.8 < z < 2.2 range 
using data from BO SS (this paper; red circles), SDSS (black squares; Richards ct al. 2006b), 2SLAQ QSO survey (light blue up-trianglcs; 
ICroom et~al . 2009a), the "boss21+MMT" survey (dark blue circles; IPalanaue-Delabrouille et al.ll2012f ) and the COMBO-17 survey (orange 
down-triangles; Wolf ct al. 2003). Center: Measurement in the 2.2 < z < 2.6 range; Right: Measu rement at z ~ 3. 2, now adding data 
from the COSMOS survey (purple diamonds; IMasters et"aT|[20T2l) and SWIRE (dark green squares; ISiana , e t al.[ |2008T). In each panel, we 
have divided the various QLFs by our best- fit "log-linear" LEDE model, which is described below in Sec. I6.ll and Table l8l 



with the transformation: 

M0 = 2) = M e (z = 2) - 0.25 (6) 
= M 1450 -0.29. (7) 

One underlying assumption in these conversions is 
that the 2SLAQ, and indeed the BOSS, quasars have 
a distribution in spectral power-law slopes (a v , where 
F(v) oc v av ) in the UV/Blue/Optical that is compara- 
ble to that of the SDSS quasar sample. Although the 
BOSS target selection avoi ds sources that would satisfy 
a UV Excess selection (see iRoss et aLll2012t iParis et al.l 
120121 ) . there is not strong a priori reason to suspect that 
these populations would deviate from a range of intrinsic 
slopes — 1 < a v < 0, centered around a„ ~ —0.40. 

Our comparison to the 2SLAQ result is shown in the 
left and center panels of Fig. Q3J for the redshift range 
1-8 < z < 2.2 and 2.2 < z < 2.6, respectively. We 
note that the 2SLAQ result is based on a combination 
of the 2SLAQ QSO survey (which dominates the signal 
at the faint end of the QLF) and the SDSS results from 
DR3 (which is responsible for the bright end measure- 
ment). Thus, the 2SLAQ points, represented as light 
blue upwards-pointing triangles in Fig. 1141 are not inde- 
pendent from the SDSS (black) squares. 

Concentrating on the z « 2.0 panel, at the bright, 
Mi(z = 2) < -26, end the boss21+MMT points, 
given by the blue filled points, seem ~0.2-0.4 dex higher 
than e.g. the 2SLAQ and SDSS data, though are gen- 
erally consistent within the quot ed (statistical) error. 
IPalana ue-Dela brouille et al.l ()2012D explore the variabil- 
ity selection in more detail than presented here, and res- 
olution of this issue will be aided by new, forthcoming, 
variability-selected data, since sample variance uncer- 
tainties over the 14 deg 2 boss21+MMT field could well 
be an issue. At the faint end, boss21+MMT, 2SLAQ and 
SDSS are all consistent. All the displayed measurements 
are consistent with the COMBO-17 points (orange down- 
triangles), due to the large error associated with those 
points (not shown). 

At z ~ 2.4 the BOSS DR9, SDSS, 2SLAQ and 



"boss2H-MMT" are in excellent agreement, at both the 
bright and faint ends. 

We compare to t he COSMOS result at z ~ 3.2 
(|Masters et al.ll2012D . in the right panel of Fig. [H The 
BOSS QLF measurement is in good agreement with the 
COSMOS results, given by the purp le upward-p o inting 
triangles. We also plot results from ISiana et al.l (|2008l 
green squares), who use an optical/infra-red selection 
over 11.7 deg 2 from t he Spitzer Wide-area Infrared Ex- 
tragalactic (SWIRE; ILonsdale et "ail 120031 ) Le gacy Sur- 
vey. The measurements from the 3 < z < 3.5 bin from 
IPalanaue- Delabro uille et all ()2012l ) is also given (blue 
circles). The BOSS DR9, SDSS, and boss21+MMT data 
are all in good agreement, and consisent given the errors. 
The faintest BOSS points are consistent with the bright- 
est SWIRE, COMBO-17 and COSMOS points, again 
given the associated errors. There seems to be an in- 
flection around Mi(z — 2) « —25.5, suggesting that our 
best-fit model is under-predicting the QLF at the both 
the bright and faint end, i.e. the bright end slope of 
the model is too steep, while the faint end slope is too 
shallow. We discuss this further in Sec. 16.11 

The R06 points lie below the other determinations, 
suggesting that they slightly unde restimated the number 
densit y of 3.0 < z < 3.5 quasars. iWorseck fe; Prochaskal 
(120111) used UV data from the GA LEX satellite 
jMartin et aLll2005t IMorrissev et al.H2007fl , to show that 
the SDSS quasar target selection systematically misses 
quasars with blue u — g < 2 colors at 3 < z < 3.5 and 
preferentially selects quasars at these redshifts with in- 
tervening Hi Lyman limit systems, causing the QLF 
to be underestimated . Inde ed, we specifically use the 
IWorseck fc Prochaskal (|2011f ) Monte Carlo model to de- 
scribe the Hi Lyman series/forest and continuum ab- 
sorption when creating our BOSS selection function, so 
we have corrected for this effect. 

6. QLF FITS, MODELS AND DISCUSSION 

In this section, we fit parametric models to our binned 
QLF and examine the evolution of the fitted parameters 
with redshift. We then compare our data to predictions 
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TABLE 7 

The binned BOSS DR9 Quasar Luminosity Function. (1) 
The mean redshift of the bin; (2) The mean s-band 
absolute magnitude of the bin; (3) the absolute magnitude 
bin center; (4) the number of quasars in each bin; (5) <£ in 
units of mpc -3 mag -1 and (6) the (poisson) error on <3>, 
divided by 1 x 10 -9 ; the bins with no measured <e> are at 
the faint end limit where the selection function is rapidly 
approaching, or is equal to, 0.00, thus making our qlf 
estimation very uncertain. however, these bins are 

INCLUDED SO THAT ^ Nq = 23 301. 



based on more physical models of quasar evolution. Fi- 
nally, we place our results in a broader context regarding 
the AGN population and its link to galaxy evolution. 

6.1. QLF model fits 

The QLF is traditionally fit by a double power-law of 
the form in Eq. [T] This functional form has four basic 
parameters, and various phenomenological models have 
been proposed to describe how those parameters evolve 
with redshift. In Pure Luminosity Evolution (PLE), 



only the break magnitude/luminosity evolves, leaving the 
overall number density constant. The opposite occurs in 
Pure Density Evolution: the shape of the QLF remains 
constant while the number density evolves. Various hy- 
brid models allow both to vary but hold the bright- and 
faint-end slopes fixed. In Luminosity Evolution and Den- 
sity Evolution (LEDE), M* (z) and $*(z) evolve inde- 
pendently, while in Luminosity Dependent Density Evo- 
lution (LDDE), the evolution of $*(z) is related to that 
of Mj*(z). Finally, extensions to these models allow the 
power law slopes to evolve as well. 

We begin with a simple PLE model for our data. In 
princi ple, M*(z) can ta ke any functional form, but we 
follow lBovle et all (|2000D by fitting it with a second order 
polynomial: 

M*(z) =M*(z = 0)- 2.5(k 1 z + k 2 z 2 ) . (8) 

We note that this quadratic form for M-* (z) requires sym- 
metric evolution about the brightest M* value, and that 
this is kn own to break down at redshifts well above the 
peak fe.g. lRichards et a l. 2006b). However, we are moti- 
vated to continue to use the quadratic PLE description as 
a historical reference and because over a limited redshift 
range the general form of our QLF is qualitatively con- 
sistent with a PLE model. For example, if the solid red 
line (representing the BOSS DR9 QLF at2.2<z<2.6) 
in Fig. [TT]is compared to the measured QLF at z < 3, 
one sees the broader trends in the data are encapsulated 
by a shift in M* with little change in normalization. 

We fit the PLE model with Eqn. |8] to our data over 
various redshift ranges. We use the combination of SDSS 
(R06), boss21+MMT and BOSS Stripe 82 dataset to per- 
form the fits. These data span 0.30 < z < 4.75 in red- 
shift, -29.55 < Mi(z = 2) < -22.96 in magnitude and 
$ = 2.2 x 10~ 9 - 2.2 x 10~ 6 Mpc~ 3 mag" 1 in number 
density. We fit to the Stripe 82 data, since we expect that 
this data is less affected by systematics, and thus more 
meaningful x 2 values can be obtained from the statistical 
uncertainties 30 . We have also found that the S82 data is 
a fair representation of the DR9 data (Fig. [T3|) . 

We perform x 2 fits to the binned data with six total 
free parameters in the PLE model, using the Levenberg- 
Marquardt optimization method to find the best-fit pa- 
rameters by minimizing the % 2 . The parameter val- 
ues for our best-fit PLE models are given in Table |S1 
We first restrict our fits to z < 2.2, where previous 
work has generally found that PLE models provide a 
reasonably good fit. Fitting over 0.30 < z < 2.20 re- 
sults in x 2 l v = 155/75. Most of the disagreement with 
our data comes at z < 1; by restricting to the range 
1.06 < z < 2.20 the fit improves to x 2 /^ = 83/52. Thus 
we find that PLE models do indeed provide a reasonable 
description of our low redshift data, though clearly there 
is room for improvement in the x 2 . At higher redshift, 
the PLE model fails. Within the BOSS redshift range 
of z = 2.20 - 3.50 we have x 2 /^ = 286/113; the re- 
sult is even worse over the full redshift range of our data 
(z = 0.30 - 3.50), with x 2 lv = 662/195. 

Fig. [TS] demonstrates why this is the case, and where 

30 Note that while Stripe 82 does not have a correction applied 
for color selection effects, the fc-correction still introduces uncer- 
tainty that may have systematic trends with redshift and luminos- 
ity. 
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the PLE model breaks down. Here we show the behavior 
with redshift of the parameters M* , a and /?. The 
parameter values and uncertainties are determined by 
y 2 minimization in each redshift bin independently. In 
the top right panel, we see that although at z < 2.2, 
a quadratic description of the evolution of M* describes 
the general trend of the data, at z > 2.2, M*{z) continues 
to get brighter and does not exhibit the turn over needed 
for the parameterization given in equation [8] to work. 
However, even if a new description of the evolution of 
M* could be found, the PLE model, with no allowance 
for density evolution, is not suitable. This is shown by 
the top left panel of Fig. [15] we can see that there is 
essentially no evolution of $* across the range 0.5 < z < 
2.2, but then log4>* declines in a roughly linear fashion 
with redshift at z > 2.2, corresponding to a drop in <E>* 
by a factor of ~6 between z = 2.2 and z = 3.5. 

Motivated by the evolution of log<I>*(z) and M*(z) 
seen in Fig. [T5] across the range z = 2.2 — 3.5, we imple- 
ment a form of the LEDE model where the normalization 
and break luminosity evolve in a log- linear manner; e.g., 

log[$*(z)] =log[$*(z = 2.2)] + Cl {z - 2.2) (9) 
M*(z)=M*(z = 2.2) + c 2 (z - 2.2). (10) 

For the BOSS Stripe 82 data across the redshift range 
z = 2.2 — 3.5, this returns a value of \ 2 l v — 136/113, in- 
dicating a reasonable fit to the data. If we instead fit to 
the BOSS DR9 data, we find generally good agreement 
in the fitted parameters, but a dramatically worse % 2 
value. While the binned QLF data from DR9 and Stripe 
82 are in good agreement, the statistical uncertainties in 
the DR9 data are far smaller due to the much greater 
number of quasars. This inflates the x 2 f° r the same 
model fit; however, as explained in Sections 13.41 and 15.11 
we expect the true uncertainties of the DR9 data to be 
dominated by systematics, in particular, in the need to 
correct for color selection effects without knowing the 
true distribution of quasar colors. The systematic effects 
associated with correcting for the selection function are 
obviously a general problem, and are especially prob- 
lematic as the selection function affects the points in a 
correlated way. Here we have taken advantage of two 
quasar samples selected by independent means; we leave 
this as a cautionary note for surveys relying on an un- 
known selection function, particularly where the data is 
dominated by objects found in regions where the selec- 
tion efficiency is low. 

We do not extend our LEDE model below z — 2.2, 
where it clearly would not describe the data. We also 
note that PDE models cannot capture the strong evolu- 
tion in M* and are easily ruled out. Qualitatively, there 
is no clear relationship between the smooth evolution of 
M* and the disjoint behavior of log$*, thus we also do 
not consider LDDE models. In summary, our data is 
best described as PLE evolution until z ~ 2.2, at which 
point a transition to LEDE evolution occurs. 

Finally, we do not see evidence for evolution in the 
power law slopes, though these are not well constrained 
by our data. In particular, we do not find a strong 
evolution of the bright end slope (bottom right panel 
of Fig. [T5"]) at z > 2.5, in contrast to R06. This could 
be because the evolution in M* affects the R06 results, 
or, very likely since we resolve the break in the QLF at 




z 

Fig. 15. — The best fit values for the parameters <3?*, M* , a and 
as a function of redshift. A double-power law model is fit to 
the Stripe 82 data in each redshift bin. The teal points are for 
the Stripe 82 data, while the red points are the BOSS DR9 CORE 
data. The four colored lines represent the four best fitting PLE 
models in Table \E\ over the respective redshift ra nges , while the 
solid black line is the log-linear LEDE model (eqn. ITOl l. 



z = 2.2 — 3.5, and consequently fit a double-power law 
model (cf. the single power- l aw in R06; see also the dis- 
cussions in lAssef et al.l[20il IShen fc Kellvl[20Ll . How- 
ever, comparing the points from Fig. 21 of R06, to our 
Fig. [15] the bright end slope measurements are consistent 
with each other, given the error bars. 

In Fig. HH we show our best-fit PLE and LEDE mod- 
els in three redshift bins, and compare our fits to other 
models that hay e been presented in the literature. 

iCrotonl (|2009l) presented a modifi cation of the PLE fit- 
ting function of lCroom et al.l (|2004f) in which the decline 
of M* with redshift is softened and the bright end power- 
law slope evolves above z = 3. This was found to fit the 
higher redshift SDSS data better than the original fitting 
form, which was fit only to the 2QZ data. We reproduce 
this modified fit in Table [5] and Fig. [TBI shown by the 
dotted (blue) line. We see that this model describes the 
data well at z ~ 2.0 and 2.4, but has a too high a nor- 
malization and (potentially) too flat a faint-end slope at 
z ~ 3.2. 

iHopkins et al.l (|2007l ) collected a large set of QLF mea- 
surements, from the rest-frame optical, soft and hard X- 
ray, and mid-IR bands, in order to obtain accurate bolo- 
metric corrections and thus determine the bolometric 
QLF in the redshift i nterval z = — 6. The observational 
dataset assembled bv IHopkins et al.l (|2007| ) is impressive, 
though most of the power in the z > 2 dataset is from 
the (R06) optical measurements of the Q LF. Taking the 
tradit ional double-power law approach, IHopkins et al.l 
(2007) then derive a scries of best-fit models to the 
QLF, including a PLE and luminosity-dependent den- 
sity evolution (LDDE) model. Their "Full" model, which 
is an LDDE-based model and includes a luminosity- 
dependent bolometric correction, is shown in Figure [12] 
by the (turquoise) dashed lines. This model fits the data 
well until the highest redshift bin at z ~ 3.2. In the 
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Fig. 16. — The BOSS DR9 Quasar Luminosity Function compared to a series of QLF fits. Left: Measurement in the 1.8 < z < 2.2 range; 
Center: Measurement in the 2.2 < z < 2.6 range and Right: Measurement at z ~ 3.2. In each panel we plot our best-fits PLE fit, given 
by the solid (purple) line, which is the fit over the redshift range 0.4 < z < 2.2 (top line in Table [8}, Also shown is our best-fit log-linear 
LEDE model, given by the (orange) dot-dashed line, with the fitting parameters al so in Table [8] The ex tension to the 2QZ QLF as given in 
ICrotonl (120091 ') . is shown by the dotted (light blue) line, while the "Full" model, of Hopkins ct al. (2007) is given by the (turquoise) dashed 
lines. 



z, redshift 
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Fig. 17. — The comparison of our best fit phenomonological mod- 
els, dashed lines, to the SDSS, BOSS Stripe 82 and boss21+MMT 
QLF data (points). The number density of various magnitude 
bins are shown as a function of cosmological time The best-fitting 
PLE model over z = 0.3 — 2.2 and best-fitting LEDE model over 
z = 2.2 — 3.5 from Table [5] are given by the dashed curves. A 
mismatch in number density at z = 2.2 for the fainter magnitudes 
is apparent, but since we do not require the fits to link, is not 
surprising and within the uncertainties. 



Hopk ins et al.l (|2007f) model, the break luminosity turns 
over at z ~ 2 and becomes fainter at higher redshift, 
while the bright end slope flattens and the normalization 
is constant. This is apparent in Fig. II 61 where the break 
luminosity is clearly much fainter than in our data and 
the faint end number densities are overpredicted. 

Fig. [17] shows the redshift evolution of the QLF in a 
series of luminosity bins, including both our data and the 
best-fit PLE+LEDE model. Previou s measurements, es- 
pecially in the deep X-ray studies (IMivaii et al.l 120011 : 
lUeda et al.ll2003t iHasinger et "ail 120051). and then in the 
optical by the 2SLA Q QSO survey ([Richards et alj|2005t 
iCroom et ail l2009al) see the trend for "AGN Downsiz- 
ing", with the number density of fainter AGN peaking 



at lower-redshift than the luminous AGN. These studies, 
especially in the optical, have generally suggested that 
PLE works up to z « 2, but not to higher z. Our BOSS 
results agree with this statement, but we use the longer 
redshift baseline of our data, and in particular the fact 
that we have resolved the break luminosity to z ~ 3.5, 
to find a simple prescription for the evolution at z > 2. 
Interestingly, we find that the shape of the QLF does 
not change (in terms of the power law slopes). Going 
from high to low redshift there is a build up of quasar 
activity (the log-linear trend in $*) until z « 2, at which 
point the number density stalls. In this LEDE-to-PLE 
toy model scenario, AGN downsizing is then simply a 
trend in L*(z). 

Our optical QLF results are also in general agreement 
with the latest determinatio n of the hard, 2-1 keV X-ray 
luminosity function (XLF; lAird et aLll2010D . These au- 
thors also find an LEDE model (which they name LADE) 
describe their XLF well, and that an XLF that also re- 
tains the same shape, but shifts in luminosity and den- 
sity, describes t he observed evolut ionary behavior. We 
also agree with lAird et al.l (|2010Q in that the (QLF) 
LEDE model shows a much weake r signature of "AGN 
Downsizing" than pre vious studies (jHasinger et alj[2005t 
iSilverman et al.ll2008D. One caveat here is that the hard 
X-ray samples used in lAird et al.l (|2010f ) are most secure 
at z < 1.2. Overall, these trends of a simple log-linear 
LEDE model describing both the QLF and XLF lends 
weight to the theory that the X-ray selected AGN pop- 
ulation at z ~ 1 is a direct descendent of the optical 
quasar population at z ~ 2; a scenario also suggested by 
quasar and X-ray AGN clustering results dHickox et al.l 
2009| IRoss et alj ' l2009l; IKoutoulidis et al.ll2012D ~ 

6.2. Quasar model predictions 

There are many models for quasar evolution in the lit- 
erature, but the modern ones come in three basic flavors. 
The first implements some of the quasar physics directly 
into numerical hydrodynamic simulations of galaxy for- 
mation or interaction. The second follows much of the 
same physics semi-analytically. The third tries to re- 
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TABLE 8 

Values from a set of the best fit double-power law evolution models (e.g. Eons. l2l l8l and ITot . Listed are the redshift 

RANGES OF THE DATA FITTED, AND THE BEST FIT VALUES OF THE MODEL PARAMETERS. We PERFORM OUR FITS USING ONLY THE 

STATISTICAL ERROR ON THE QLF. 
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late the properties of quasars and black holes directly to 
those of dark matter halos or the galaxies which reside in 
them. We give recent examples from each of these classes 
of models here. 

DeGraf et al. (2013, in prep.) present models for the 
QLF using the new "MassiveBlackH" hydrodynamic sim- 
ulation, which has a boxsize of 100 hr 1 Mpc, number 
of particles, iV p = 2 x 1792 3 and a gravitational soft- 
ening of e = 1.85 h " 1 kpc, and employs a WMAP7 
IjKomatsu et al.lfeOllD cosmology. These simulations in- 
corporate the physics of hydrodynamics, radiative cool- 
ing, star formation, black holes and associated feedback 
in order to make ab initio predictions for the observed 
properties of galaxies and quasars. The QLF for each 
redshift bin is computed using the complete luminosity 
history of every black hole, producing the best available 
statistics and extending the predictions to the brightest 
luminosity by catching rare objects that only occasion- 
ally reach very high-L. The predictions from these hy- 
drodynamic simulations are given by the shaded black 
region in Fig. [T5J Note that they extend to luminosities 
fainter than BOSS generally probes. 

There are discrepancies between the simulations and 
the data, especially at z ~ 2.0 and 2.4, which may be 
due to severa l effects. Previous work on smaller sim- 
ulations (e.g.. iDegraf et aL"ll2010D found that lower res- 
olution simulations produce steeper faint end luminos- 
ity functions. Thus increased resolution should further 
flatten the faint end. At the bright end, volume limi- 
tations become significant, with only several black holes 
reaching the brightest luminosities. The shaded region 
in the bright end QLF represents an estimate for the cos- 
mic varianc e using a larger volume simulation ( "Massive - 
Black" , see IDeGraf et all [2011 IDi Matteo etaTl [20121 ) . 
The larger simulation avoids the volume limitations re- 
sulting in the upper bound of this region, suggesting that 
within volume limitations the simulations are consistent 
w ith current da ta. 

iMarulli et al.l (|2008| ) model the cosmological co- 
evolution of galaxies and their central supermassive black 
holes within a semi-analytical framewor k developed on 
the o utputs of the Millennium Simulation (Springcl et al. 
20051). These autho rs use the galaxy formation model of 
Croton et al.1 (|2006l ) as updated by iDe Lucia fc Blaizotl 
( 20071 ) as their starting point. Luminous quasars in this 
model occur when a BH accretes cold gas after a major 



merger of two gas rich galaxies. The accreted mass is 
proportional to the total cold gas mass present, but with 
an efficiency which is a function of the size of the sys- 
tem and the merger mass ratio, and chose n to reproduce 
the ob served local Mbh — Mb u i ge relation. IMarulli et al.l 
(2008) then couple this accretion to various light curve 
models. The predictions for the luminosity function are 
shown in Fig. [18] by the triple-dotted-dashed (turquoise) 
line. We see that this model does well in the lower red- 
shift bins at z ~ 2.0 and 2.4 at reproducing the data, 
but perhaps over predicts the number of faint quasars at 
z ~ 3.2. 

For c omparison we also cons ider a second semi-analytic 
model (|Fanidakis et al.l [20121 ). This model is embed- 
ded in t he semi-analytica l galaxy f ormation cod e GAL - 
FORM (ICole et aL[ l2000l see also IBaugh et all (pOOl ; 
IBower et alj (120061)) and predicts the masses, spins 
(|Fanidakis et al.l 1201 If) and mass accretion histories of 
BHs in tandem with the formation of their host galax- 
ies. In addition to merger-induced triggering they al- 
low triggering when discs bec oming dynamically un stable 
(b ased on the a r gume nts in Efstathi ou et al.l I1982T) . As 
in lMarulli et al.1 (|2008[ ) they also follow quasi-hydrostatic 
hot gas accretion (known variously as "hot halo mode" , 
"radio mode" or "radiative mode" accretion) with a rate 
orders of magnitu de below the Eddingt on limit. The 
key aspect of the Fanidaki s et al.l ([2012) model in our 
comparison is that their starburst mode, and thus the 
BH mass grow th, is mainly driven by disc instabilities. 
Comp arison of IMarulli et al.l (|2008|) and iFanidakis et all 
( 2012) thus allows insights into how the triggering mode 
of quasar activity can potentially be tested by measure- 
ments such as o urs. The number densities from the 
IFanidakis et al.l ()2012l ) model are calculated consider- 
ing the entire population of AGN (both obscured and 
unobscured) a nd includ e the e mpirical obscuration pre- 
scription from [Hasinger (200§j). The QLFs for the unob- 
scured population are shown as (purple) dot-dashed lines 
in Fig. [H 

Hirsch mann et all (I201 2 J) also used semi-analyt ic mod- 
els, based on those from ISomerville to ex- 
amine the properties of accreting B Hs and the evolu- 
tion o f the QLF. (We do not show the Hirschmann et al. 
( 2012) predictions in Fig. [181 but their best fitting model 
fits our data well with potentially a slight overproduction 
of the faintest QSOs at z > 2.5; see their Fig. 7.) These 
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Fig. 18. — The BOS S DR9 Quasar Luminosity Function compared to a series of QLF models from the literature. The model from 
Conro^ & White (2012j) is given by the dashed orange line, while th e model fro m Shcn (2009) is given by the light-blue dotted line. The 
MaruIIi et al . (2008) model is given by the solid turquoise line, the Fanidakis ct al. (2012) model by the dot-dashed purple line and the 
model from (DeGraf et al. 2013) is the shaded black region. We refer the interested reader to the given papers for presentation and 
discussion of the uncertainties associated with the published models. Left: Measurement in the 1.8 < z < 2.2 range; Center: Measurement 
in the 2.2 < z < 2.6 range; Right: Measurement at z ~ 3.2. Note, the Shen09 and CW12 models are on top of each other at Mi < —26. 



authors find that their best fitting model (which includes 
using "heavy" black hole seeds of M soo d ~ 1O 5_6 M0 
at very high z and a varying sw&-Eddington limit for 
the maximum accretion rate at z < 1) suggests a sce- 
nario in which the disc instabilities are the main driver 
for moderately luminous Seyfert galaxies at low redshift, 
but major mergers remain the key trigger for luminous 
A GN/quasars , especially at high z. 

Shcn (2009) presents a phenomenological model for 
the growth and cosmic evolution of SMBHs, in which 
the quasar properties are tied to the properties of dark 
matter halos, rather than galaxies drawn from a semi- 
analytic model. This model assumes that quasar activity 
is triggered by major mergers of host halos, and that the 
resulting light curve follows a universal form, in which 
its peak luminosity is correlated with the (post)merger 
halo mass. Quasar activity is quenched at low z and in 
lower mass halos with phenomenological rules. In partic- 
ular, the quasar trigge ring rate de pends on a "quasar-on" 
factor (called /qso m lShen| [2009') which has exponential 
cut-offs both at the low and high mass ends which are 
adjusted to fit the data. These cut-offs ensure that halos 
with too small a (postmerger) halo mass cannot trigger 
any quasar activity, while those above a (redshift depen- 
dent) maximum mass cannot cool gas efficiently and BH 
growth halts. With these assumptions, the quasar LF 
and SMBH growth are tracked self-consistently across 
cosmic time. The QLF predicted by this model is shown 
in Fig. [TS] by the dotted (blue) line. This model does 
well at reproducing the data in all three redshift slices, 
though with a slight over-production of bright quasars at 
z~2.0. 

Recently IConrov fc White! (|2012D presented a model 
for quasar demographics in which quasars populate 
galaxies in a simple manner and many of the proper- 
ties of the quasar population follow naturally from the 
known, evolving properties of galaxies. A simple "scat- 
tered lightbulb" model is adopted, with BHs shining at a 
fixed fraction of the Eddington luminosity during accre- 
tion episodes with Eddington ratios drawn from a log- 



normal distribution. The quasar duty cycle is explicitly 
independent of galaxy and BH mass and l umino s ity, in 
contrast to the strong dependence invoked in lShenl (2009) 
when connecting quasars to halos. The QLF predictions 
for that model are shown in Fig. [18] as the (red) dashed 
lines. 

While the models we have highlighted agree with the 
existing data relatively well, they explain the qualitative 
behaviors we see in different ways. For example, it is 
well known that the abundance of bright quasars drops 
rapidly to low z and that lower mass black hole growth 
peaks at lower redshift than higher ma ss black holes 
(IHasinger et al . 2005; C ropm et a l. 2009a). In the model 
of IConrov fc White] (|2012f ) this is explained through a 
combination of slow growth of massive galaxies an d evo- 
lution in the Eddington ratio. In the model of iShenl 
(2009), it involves a suppressing function which simu- 
lates the effec ts of cold gas consump tion with time. In 
the model of iFanidakis et all ([2012) it arises due to a 
combination of factors, including obscuration evolution. 

The models differ significantly in the mass and red- 
shift dependence of the duty cycle, and predict subtle 
differences in the width of the halo mass distribution 
at any redshift. In almost all models the characteristic 
halo mass associated with existing quasar samples is al- 
most independent of redshift. This arises largely due to a 
chance cancellation of trends in the absolute magnitude 
limit, the relation between galaxy and halo properties 
and galaxies, blac k holes and Eddington ratios. 

In the model of IConrov fc White! (|2012fl , the evolution 
of the characteristic luminosity is driven by the evolu- 
tion in the L — M ga i and Af ga i — Mh relations, while the 
break in the LF arises pr imarily due to the shape of the 
Mgai - M h relation. The IShenl ()2009D model adjusts the 
typical host halo of luminous quasars to fit the observed 
evolution of the break luminosity. In the semi-analytic 
models, the starburst/quasar mode is powered by one 
or a combination of major galaxy mergers and disk in- 
stabilities with the relative contributions possibly evolv- 
ing with time. The evolution of the characteristic lumi- 
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nosity thus arise s from a complex interplay of factors. 
iConrov k White] ()2012[ ) predict that the faint end slope 
of the LF does not vary significantly, while the bright 
end s l ope ap pears shallower at higher z. In the model of 
Shen (2009), the LF is predicted to turn down at suffi- 
ciently low luminosities and high redshifts, since at z > 2 
the minimum Eddington ratio is constrained by the age 
of the Universe. The hydrodynamic simulations predict 
a steep faint-end slope at z ~ 2. Most of the models 
have considerable scatter between quasar luminosity and 
galaxy or halo mass, and thus predict a power-law tail 
to high luminosity, as observed. Further measurements 
of this tail at higher z may provide better constraints on 
this aspect of the models. 

In summary, all the models reproduce the QLF and 
quas ar demographics overall reasonably well. We agree 
with iHirschmann et all (|2012| ) when they state that fur- 
ther progress on these issues will require data beyond 
just the luminosity function. 

6.3. Discussion 

In this final section, we tie our QLF results (and com- 
parisons to models) into the broader context of the link 
between SMBH growth (see as seen AGN activity) and 
the properties of galaxies. We take as our starting point 
the QLF reported here and the clustering measurement 
and discussion of the BOSS DR9 uniform quasar sam- 
ple reported i n I White et all (120121) . Using the same ar- 
gu ments as in iWhite et aLT(l2012|) 1 and the c onversions 
of lCroom et all (|2005| ) and lShen et al.1 (|2009l ). we place 
the median BOSS quasar with a bolometric luminosity of 
Lbol = 2 — 4 x 10 46 erg s _1 , in dark matter haloes of char- 
acteristic mass of ~ 2 x 10 12 h~ x M@ at z « 2.5. Either 
making the assumption that the BOSS q uasars are con- 
sistent with the Mbh — Mh a i relation (jFerraresel 120021 : 
iFine et all 120061 ). or, that the quasars radiate at close 
to the Eddington Limit, L Edd = 10 47 1 (M B h/10 9 M q ) 
erg s" 1 , suggests that the median Mbh in our sample 
is - 2 x 10 8 M Q . As a guide, a typical - 2 x 10 8 M Q 
BH, accreting continuously since z ~ 2.5, with an ac- 
cretion efficiency of e = 0.1, and not merging, would 
have a mass at redshift z ~ of Mbh ~ 6 x 1O 1O M . 
This would place these objects at the very highest BH 
masses observed, but also inline with recent results 
(|McConnell et alJ20ll . A more realistic scenario, where 
the duty cycle is 1%, would lead to Mbh ~ 6 x 1O 8 M , 
placing these objects in bulges with a ~ 250 — 300 km 
s -1 , and thus in early-typ e galaxies from the relations in 
e.g. iGultekin et ail (|2009l ). 

From the observed clustering (and indeed essentially 
any of the models quoted above) the typical halo for a 
BOSS quasar at z ~ 2.5 would grow to host a small 
group by z ~ 0. The most likely host galaxy is the cen- 
tral galaxy of the group, since at higher- z, any satellites 
would not be massive enough to host a SMBH. Thus, 
the typical BOSS quasar host descendant would be the 
central galaxy of a small group - though we caution that 
including e.g. the diversity of growth histories of DM 
halos and scatter in any of the given relations, can easily 
lead to an order of magnitud e dispersion in the above 
statements ()White et alll2012l ). Placing these quasars at 
the centers of groups at z = is consistent with the sug- 
gested velocity dispersions given above. This potentially 



also suggests that BOSS quasars today are very likely 
not on the "SF Main Sequence" any more (i.e. they are 
quenched) even if they were i nitially. This i s also consis- 
tent with the recent work bv lKellv k Shenl (|2012t ). 

Leaving the properties of the median BOSS quasar, we 
now focus on the "extremes" of our population. Taking 
the most luminous quasars, we find these objects to have 
close to log(-Lboi) = 46.0, and thus black holes in the mass 
range ~ 3 x 10 9 M Q (assuming an Eddington luminosity). 
At the bright end, the QLF is described by a power-law 
fall-off, while the massive end of the stellar mass func- 
tion, the abundance declines exponentially. With a rela- 
tionship known to exist between M-BB./M\,ulge (and where 
Mbuigo ~ Mg a i for these compact massive galaxies), this 
argues that there is scatter in Lq at fixed M m \. This is 
perhaps not surprising: at low z, MsH/Mbuige is mea- 
sured to have ~0.3 dex in scatter and Eddington ratios 
are also measured to have ~0.3 dex scatter, so a scatter 
of at least 0.4 dex overall could be expected. However, 
this leads to the situation that at high-L, scatter is in- 
creasingly important, and that bright quasars are "over- 
bright" , and it is currently unclear what underlying phys- 
ical mechanisms would lead to this enhanced up-scatter. 
We leave further investigation into the potential evolu- 
tion of MBH/M ga i, and the different channels that drives 
the growth of black holes, the evolution of the number 
density of quasars, and that of AGN activity in general 
for future study. 

7. CONCLUSIONS 

The quasar luminosity function is one of the most fun- 
damental observables of this class of important cosmolog- 
ical objects. The shape and evolution of the QLF pro- 
vides constraints on models of quasar fueling, feedback 
and galaxy evolution and the ionization history of the 
inter-galactic gas. Despite its importance, it has proven 
difficult observationally to probe the quasar luminosity 
function at magnitudes below the break at the peak of 
the quasar epoch. 

Here we measure the QLF using data from the SDSS- 
III: Baryon Oscillation Spectroscopic Survey (BOSS) us- 
ing a uniformly selected sample of 23 301 quasars, and fill 
in the L — z plane with published results from the SDSS- 
I/II. We probe the faint end of the QLF to M l = -24.5 
at z — 2.2 and complement our uniform color-selection 
with a sample of variability-selected quasars from the 
"Stripe 82" field. We also provide a cross-check of our 
selection function using new, simulated, model, quasar 
spectra. Amongst our findings are: 

• That down to a magnitude limit of i = 21.5, there 
are 26.2 and 48.0 quasars deg -2 across the redshift 
ranges 2.2 < z < 3.5 and 1.0 < z < 2.2 respec- 
tively. Using the deeper boss21+MMT data, for 
the unobscured 1.0 < z < 2.2 quasar population, 
there are 78 objects deg" 2 brighter than i sa 23.0, a 
surface density similar to that selected by a shallow 
mid-infrared selection (|Stern et a l. 2012). 

• Our combined SDSS+BOSS QLF is reasonably 
well described by a double power-law, quadratic, 
pure luminosity evolution (PLE) model across the 
redshift range 0.3 < z < 2.2, with a bright end 
slope -3.37lg$j, a faint end slope — 1.16±g:g|, 
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M*(z = 0) = -22.8518:5$, h = 1.241°;^, fc 2 = 
-0.251^2 and log$* = -5.96l°;°g. 

• The simple PLE model breaks down at z > 2.2. 
We replace it with a luminosity evolution and den- 
sity evolution (LEDE) model that has a log-linear 
trend in both <£>* and L* . This simple form pro- 
vides a good fit to the data at 2.2 < z < 3.5, 
capturing both the steep decline in number den- 
sity and the rise in the break luminosity. The data 
are consistent with no evolution in the power law 
slopes, though do not strongly constrain the lack 
of evolution. 

• We compare our measured QLF to theoretical mod- 
els and find a wide variety of models describe our 
data reasonably well. While the latest hydrody- 
namic simulations do not fit as well, semi-analytic 
models in which luminous quasar activity is trig- 
gered by major mergers, disk instabilities or a com- 
bination of channels can fit our data over a wide 
range of redshifts. Models based on directly popu- 
lating halos with quasars can fit the shape of our 
QLF by assuming a mass and redshift-dependent 
duty-cycle which is sharply peaked around a char- 
acteristic mass. We also find that models which 
relate black hole mass linearly to galaxy mass and 
assume a mass-independent duty-cycle match our 
QLF well. 

The results presented here are from the first two, of 
five, years of BOSS spectroscopy. The upcoming Data 
Release Ten dataset will cover ~ 7000 deg 2 , include 
~ 150, 000 quasars and will more than double the num- 
ber in our uniform selection. Future investigations will 
be able to use this enhanced dataset in order to fur- 
ther quantify, and refine, the selection function for the 
2.2 < z < 3.5 quasar sample and thus reduce the errors 
further. This release will include quasars that were ob- 
served by BOSS because of their near- and mid-infrared 
colors, and with these samples we will be able to infer 



further key properties of quasars at the height of the 
quasar epoch. 

The JavaScript Cos mology Calcu lator was used whilst 
preparing this paper ([Wright 2006) . This research made 
use of the NASA Astrophysics Data System. Heavy use 
was made of the RED IDL cosmolog y routines w ritten 
by L. and J. Moustakas, and based on [Hogg (1991). 
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APPENDIX 

APPENDIX A. COMPARISON OF QUASAR SPECTRAL MODELS 

In Section 13.41 we introduced three models for quasar spectral features. In brief, the fiducial model (adopted for 
our primary analysis) includes a luminosity-dependent emission line template derived from fitting of composite quasar 
spectra. The composite spectra are created within narrow bins of luminosity, so that mean trends of emission line 
features with luminosity are repr oduced; in par ticular, the anti-correlation of line equivalent width with continuum 
luminosity (the "Baldwin Effect" . lBaldwiniri977l) . Introducing this feature accounts for the luminosity dependence of 
quasar colors and the effect this has on quasar selection. 

For comparison, we include two additional models. The first is based on a fixed emissi on line template with 
no luminosity depe ndence. As this is most similar to models used in previous work (e.g., iRichards et al.l l2006bl : 
iCroom et al.l l2009af ). it provides a reference point for comparison to QLF estimates that did not include a Baldwin 
Effect in the selection function estimation. Finally, we also include a model with dust extinction ("exp dust") , 
motivated by observations of SDSS quasars with mild dust reddening (e.g.. IRichards et al"1l2003t iHopkins et al.|[2004l) . 

In Figure[T2lwe compare the estimated QLFs derived from each of the three selection function models. The systematic 
effects resulting from imperfect knowledge of the true selection function (or by corollary, the intrinsic distribution of 
quasar spectral features) is greater than the statistical uncertainties resulting from Poisson variations. 

Here, we motivate our choice of selection function (the fiducial model). Our method is a simple qualitative comparison 
of the observed color-redshift relation for the three models as compared to the data. The method for constructing the 
color redshift relation is described in Section 13.41 

The resulting relations are shown in Figure [19j The "VdB lines" model does poorly at reproducing the observed 
colors in the range 2.4 < z < 3.3, when Lya and C IV are in the g and r bands, respectively. The only difference 
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TABLE 9 

The narrowly binned BOSS DR9 Quasar Luminosity Function. The columns are the same as Table [T] The full table 

appears IN the electronic edition OF The Astrophysical Journal. 



between this model and the fiducial model is the emission line template, thus a model that does not account for the 
Baldwin Effect will have difficulty r eproducing quasar col ors at these redshifts. Note that this effect is likely less 
pronounced in the SDSS data (e.g., iRichards et al.1 l2006b[ ) . as it covers less dynamic range in luminosity, and thus 
most quasars are closer to the mean luminosity represented by a single composite spectrum. The exp dust model 
appears to do as well as the fiducial model. We chose not to include dust in order to remain more consistent with 
previous work, and since the focus of this work is on the unobscured quasar population. Unsurprisingly, the exp dust 
model results in a lower overall completeness, so that the estimated luminosity function is higher overall (Figure IT2"T) . 
We will consider these issues further in subsequent work. 
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Fig. 19. — Color-redshift relations for the three quasar spectral models. Each column of panels represents one of the quasar models 
named above the top panel. As in Fig. [5] the solid blue lines are the mean and ±ltr scatter of the colors in redshift bins of Az = 0.05. The 
dashed red lines are the same color relations derived from the simulated quasars. 



APPENDIX B. ADDITIONAL BOSS QLF TABLES 



Here we present additional tables reporting the BOSS QLF for the various samples given in the main text. Tabled 
gives the BOSS DR9 QLF as shown in Fig. [TSl while Table [TOl gives the calculated QLF for the Stripe 82 dataset over 
the same redshift range and binning (teal points in Fig. IT3|) . 
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TABLE 10 

The narrowly binned BOSS Quasar Luminosity F unct ion using data from 5731 (5476) 2.20 < z < 4.00 (3.50) quasars selected 

VIA THEIR VARIABILITY SIGNATURE ON STRIPE 82 CSeC. 12.4) . THE COLUMNS ARE THE SAME AS TABLE THE FULL TABLE APPEARS IN THE 

ELECTRONIC EDITION OF The Astrophysical Journal. 
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